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CN ! Abstract 

q ■ In this paper we present a new characterization of Sobolev spaces on W\ Our 

characterizing condition is obtained via a quadratic multiscale expression which ex- 
ploits the particular symmetry properties of Euclidean space. An interesting feature 
■ of our condition is that depends only on the metric of R" and the Lebesgue measure, 

so that one can define Sobolev spaces of any order of smoothness on any metric 
^ , measure space. 

u 

"5 : 1 Introduction 

In this paper we present a new characterization of the Sobolev spaces W a ' p on R", where 
the smoothness index a is any positive real number and 1 < p < oo. Thus W a ' p consists 
of those functions / e U = L p (R n ) such that (-A) a/2 f e LP. Here A is the Laplacean 
(>. \ and (—A) a/ f is defined on the Fourier transform side by \^\ a f{^). If < a < n this means 

that / is a function in U which is the Riesz potential of order a of some other function g 
in L p , namely / = c n l/\x\"~ a * g. If a is integer, then W a,p is the usual space of those 
functions in U such that all distributional derivatives up to order a are in LP . 

To convey a feeling about the nature of our condition we first discuss the case a = 1. 
Consider the square function 



fB( X ,t)- fix) 1 dt 



j, xe W . (1) 



Here / is a locally integrable function on R" and fB( X ,t) denotes the mean of / on the open 
ball with center x and radius t. One should think of ^ BU,) f{x) as a quotient of increments 
of / at the point x. Our characterization of W l ' p reads as follows. 



Theorem 1. Ifl < p < oo, then the following are equivalent. 

(1) f€ W hp 

(2) feL p andS(f)eL p . 

If any of the above conditions holds then 

lis(/)ll p -liv/v 
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The symbol A B means, as usual, that for some constant C independent of the 
relevant parameters attached to the quantities A and B we have C l B<A<CB. 

Notice that condition (2) in Theorem Q] above is of a metric measure space character, 
because only involves integrals over balls. It can be used to define in any metric measure 
space X a notion of Sobolev space W l,p (X). It is not clear to the authors what are the 
relations of this space with other known notions of Sobolev space in a metric measure 
space, in particular with those of Hajlasz flHll or Shanmugalingam (Si (see also HHKIO . 

The proof of Theorem Q] follows a classical route (see UStrlO . The relevant issue is the 
necessary condition. First, via a Fourier transform estimate we show that 



for good functions /. In a second step, we set up a singular integral operator T with values 
in L 2 (dt/t) such that 



The kernel of T turns out to satisfy Hormander's condition, so that we can appeal to a well 
known result of Benedek, Calderon and Panzone HGRl Theorem 3.4, p. 492] on vector 
valued Calderon-Zygmund Theory to conclude the proof. The major technical difficulty 
occurs in checking Hormander's condition. 

The proof extends without pain to cover orders of smoothness a with < a < 2. The 
square function S (/) has to replaced by 



The result is then that, for < a < 2, / e W a ' p is equivalent to / e If and S a (f) e LP. 
Notice that 



where the barred integral on a set stands for the mean over that set. Stricharzt ( UStrl ) used 
long ago the above square function for < a < 1 to characterize W a,p . However the 
emphasis in UStrM was on a larger variant of S a (f) in which the absolute value is inside the 
integral in y in ©. In the interval 1 < a < 2 putting the absolute value inside the integral 
destroys the characterization, because then one gives up the symmetry properties of R". 
For instance, S a (j) vanishes if / is a first degree polynomial. 

There are in the literature square functions very close to © which characterize W a ' p , 
for < a < 2. For example, first differences of / may be replaced by second differences 
and the absolute value may be placed inside the integral ( UStrM and (Stl Chapter V]). The 
drawback with second differences is that they do not make sense in the setting of metric 
measure spaces. See also the paper by Dorronsoro llDll . 

We now proceed to explain the idea for the characterization of W 2 ' p . Take a smooth 
function / and consider its Taylor expansion up to order 2 around x 



P(/)|| 2 = c||V/|| 2 , 



\\T(f)\\i?®«,m m = \\s<j)h- 





(2) 



f(y) = fix) + V/(x) • (y - x) + ^ dPf(x)(y -xf + R, 



(3) 



iei=2 
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where R is the remainder and B a multi-index of length 2. Our goal is to devise a square 
function which plays the role of S i(/) (see © for a = 1) with respect to second order 
derivatives. The first remark is that the mean on B(x, t) of the homogeneous polynomial 
of degree 1 in © is zero. Now, the homogeneous Taylor polinomial of degree 2 can be 
written as 

da fix) a It 

^-(y -xf = H(y -x) + —Af(x) \y - x\ 2 , (4) 
B\ 2n 



for a harmonic homogeneous polynomial H of degree 2. Hence the mean on B(x, t) of the 
homogeneous Taylor polinomial of degree 2 is 



j 

JB(x, t) 



^Af(x) \y-x\ 2 dy. 
2n 



This suggests defining 



Siimy 



= r f 

Jo Jb(xJ) 



(fiy)-fix)-±iAf) B(XJ) \y-x\ 2 ) 



dy 



dt 
7' 



x e 



(5) 



We cannot replace (A/) B(t J} by Af(x) in the preceding definition, because the mean guar- 
antees a little extra smoothness which one needs in a certain Fourier transform computa- 
tion. Notice that, according to the remarks made before on the mean on the ball B(x, t) 
of the homogeneous Taylor polynomials of degrees 1 and 2, in the expression above for 
Si{f){x) one may add the missing terms to get the full Taylor polynomial of degree 2, 
except for the fact that Af(x) should be replaced by (A/) B ( V ,). Were / smooth enough, 
one could even add the homogeneous Taylor polynomial of degree 3, because it is odd 
(taking x as the origin) and thus its mean on B(x, t) vanishes. This explains why whatever 
we can prove for a = 2 will also extend to the range 2 < a < 4 by defining 



s a (f)(xy 



= r f 

Jo JBi 



(f(y)-f{x)-± n {Af) B{xJ) \y-x\ 2 ) 



dy 



>B(x,t) 



dt 

7' 



X 6 



(6) 



Here is our second order theorem. 



Theorem 2. If I < p < oo, then the following are equivalent. 

(1) f G W 2 >P 

(2) f e U and there exists a function g e U such that S2CA g) s L p , where the square 
function S 2(f, g) is defined by 



S 2 (f,g)(x) 2 



r f 

Jo Jb 



(f(y) - fix) - gB{ X ,t) \y - M 2 ) 



dy 



'B(xJ) 



dt 

— , x e 



Iffe W 2 ' p then one can take g = Af/2n and if (2) holds then necessarily g = Af/2n, 



a. e. 
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If any of the above conditions holds then 



ws{f,g)\\ p ^\m\ p . 

Notice that condition (2) in Theorem [2] only involves the Euclidean distance on W and 
integrals with respect to Lebesgue measure. Thus one may define a notion of W 2,p (X) on 
any metric measure space X. For more comments on that see section 4. 

Again the special symmetry properties of R" play a key role. For instance, S2 anni- 
hilates second order polynomials. Theorem |2] has a natural counterpart for smoothness 
indexes a satisfying 2 < a < 4. The result states that a function / G W a,p if and only if 
/ 6 LP and there exists a function g e LP such that S a (f, g) e LP, where 



s a (f,g)(xY 



r f 

JO JB(x, 1 



(f(y) - f(x) - gB( X ,t) \y - M 2 ) 



dy 



dt 

7' 



x e 



We proceed now to state our main result, which covers all orders of smoothness and 
all p with 1 < p < 00. Before it is convenient to discuss the analogue of © for homoge- 
neous polynomials of any even degree. Let P be a homogeneous polynomial of degree 2 j. 
Then P can be written as 

P(x) = H(x) + A j P — \x \ 2J , 

where Lj = A j (\x \ 2j ) and H satisfies A j H = 0. This follows readily from (Sfl 3.1.2, p. 69]. 
Considering the spherical harmonics expansion of P(x) we see that f , H{x) dcr = 0, 

J\x\— 1 

cr being the surface measure on the unit sphere, and thus that f H(x) dx = 0, t > 0. The 
precise value of Lj, which can be computed easily, will not be needed. 

Theorem 3. Given a > choose an integer N such that 2N < a < 2N + 2. If 1 < p < oo, 
then the following are equivalent. 

(1) f€ W a ' p 

(2) f g L p and there exist functions g j G LP, 1 < j < N such that S a {f,g\,g2, ■■■ ,gN) £ 
LP , where the square function S a (f, g\,gi,..., gN) is defined by 



S a (f,gl,g2,---,gN)(x) 



Jo JB(x,t) 



R N (y, x) 



dy 



dt 

7' 



X G 



and RNiy, x) is 



Rn(y, x) = f(y) - f(x) - gi(x) \y - x\ z + ■ ■ ■ + g N -i(x) \y - x\ 2(N r) + (g N ) B ( x , t )\y - x 

If f G W a ' p then one can take gj = A j f/Lj and if (2) holds then necessarily gj 
A j f/Lj, a. e. 

If any of the above conditions holds then 



.\2N 



\\S a (f,gi,...,g N )\l 



i-&T /2 f\\ P 
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Again condition (2) in Theorem [2] only involves the Euclidean distance on R" and 
integrals with respect to Lebesgue measure. Thus one may define a notion of W a ' p (X) for 
any positive a and any 1 < p < oo on any metric measure space X. For previous notions 
of higher order Sobolev spaces on metric measure spaces see ULLWL See section 4 for 
more on that. 

The proof of Theorem [3] proceeds along the lines sketched before for a = 1. First we 
use a Fourier transform computation to obtain the relation 

\\S a (f, Af/L u A N f/L N )\\i = c ||(-Ar /2 /ll2 . 

Then we introduce a singular integral operator with values in L 2 (dt/t 2a+i ) and we check 
that its kernel satisfies Hormander's condition. 

The paper is organized as follows. In sections 1, 2 and 3 we prove respectively The- 
orems [TJ |2] and [3] In this way readers interested only in first order Sobolev spaces may 
concentrate in section 1 . Those readers interested in the main idea about jumping to orders 
of smoothness 2 and higher may read section 2. Section 3 is reserved to those interested in 
the full result. In any case the technical details for the proof of Theorem \T\ are somehow 
different of those for orders of smoothness 2 and higher. The reason is that Horman- 
der's condition involves essentially taking one derivative of the kernel and is precisely the 
kernel associated to the first order of smoothness that has minimal differentiability. 

Our notation and terminology are standard. For instance, we shall adopt the usual 
convention of denoting by C a constant independent of the relevant variables under con- 
sideration and not necessarily the same at each occurrence. 

If / has derivatives of order M for some non-negative integer M, then V M f = (cPf)\p\ = M 
is the vector with components the partial derivatives of order M of / and |V M /| its Eu- 
clidean norm. 

The Zygmund class on W consists of those continuous functions / such that, for some 
constant C, 

I/O + h) + f(x -h)- 2/0)| < C \h\, x, h e R" . 

The basic example of a function in the Zygmund class which is not Lipschitz is f(x) = 
\x\ log|x|, x e R". 

The Scharwtz class consists of those infinitely differentiable functions on R" whose 
partial derivatives of any order decrease faster than any polynomial at oo. 



2 Proof of Theorem U 

The difficult part is the necessity of condition (2) and we start with this. 
As a first step we show that 

Pi(/)|| 2 = c||V/|| 2 (7) 

for a dimensional constant c. Set 

x(x) = \mv\ xmi){x) 
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and 

1 x 
Xt(x) = -Xi~) , 

so that, by Plancherel, 

f S l (f)(x) 2 dx = f f \(f* Xt ){x)-f{x)\ 2 dx^ 
Jr" Jo Jr" ' 

= c £ £m{)-M 2 \m\ 2 dfj. 

Since x ls radial, x(%) = ^(1^1) for a certain function F defined on [0, oo). Exchange 
the integration in di; and dt in the last integral above and make the change of variables 
r = ?|£|.Then 



i(F(r)-n 2 ^ \m\ 2 \t\ 2 dt 



= c f |(F(r)-l| 2 % IIV/II 2 
Jo T 



and © is reduced to showing that 



|(F(r)-l| 2 4< TO - 

o 



r 3 



Set B = 5(0, 1) and e x = (1,0, ... ,0) e R". Then 

F(t) =x(tei) = T exp (-ix\t) dx 
Jb 



ix\t - ^^i^ 2 + " " " 



= 1 - - \ x\ dx t 2 + • • • , 



which yields 

Fit) - 1 = 0(t 2 ), as t -> 

and shows the convergence of © at 0. 

Since F{\^\) = x(0 is the Fourier transform of an integrable function, F(r) is a 
bounded function and so the integral ([8]) is clearly convergent at oo. 

We are left with the case of a general p between 1 and oo. If / e W l,p , then / = 
g* 1 /\x\ n ~ l for some g e LP (with 1 l\x\ n ~ l replaced by log |x| for n = 1). Set /(x) = 1 . 
Then 

fm*i> ~ fix) = (f * Xt ){x) - fix) = (g * K t )(x) , 

where 

K t (x) = (I * xt)(x) - Kx) = f I(y) dy - I(x) . (9) 

JB(x, t) 



6 



If we let T(g)(x) = (g * K t )(x), x e K", then one can rewrite S i(f)(x) as 



S i(f)(x) = 11 \(g* K t )(xf ^ I = \\Tg{x)\\ mdtli . 



3 -II- 1 £WIIL 2 (A/f 3 ) • 



Then ([7]) translates into 



f \\Tg(x)\\ 2 LHdt/t3) dx = c\\g\\l, 
Jw 



and we conclude that T is an operator mapping isometrically L 2 (R") into L 2 (R", L 2 (dt/ 1 3 )). 
If the kernel K,(x) of T satisfies Hormander's condition 

f - y) - K t (x)\\ L 2 {dm <C, ye R" 

then a well known result of Benedek, Calderon and Panzone on vector valued singular 



integrals (see [GR, Theorem 3.4, p. 492]) yields the LP estimate 



I 



\\ T Six)\\" dx<C p \\g%, 



'L 2 (dt/t 3 ) 

which can be rewritten as 

I|Si(/)IIp<C p ||V/|| p . 

The reverse inequality follows from polarization from ([7]) by a well known argument 
( PGRi p. 507]) and so the proof of the necessary condition is complete. We are going to 
prove the following stronger version of Hormander's condition 

\\K t (x-y)-K t (x)\y idm <C-^, veR", (10) 

for almost all x satisfying \x\ > 2\y\. 

To prove (TTOl) we deal separately with three intervals in the variable t. 

Interval 1: t < j. From the definition of K t in © we obtain 

VK t {x) = (VI* X t)(x)-VI(x). (11) 

Notice that, in the distributions sense, the gradient of / is a constant times the vector 
valued Riesz transform, namely 

V/ = -in - l)p.v. 



If \x\ > 2\y\, then the segment [x - y, x] does not intersect the ball 5(0, \x |/2) and thus 

\K t (x-y)-K t (x)\<\y\ sup |V^ f (z)[ • (12) 

zs[x-y,y] 
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If t < |x|/3 and z e [x - y,y], then B(z, t) c R" \ 5(0, \x |/6), and hence 

VK t (z) = f (V/(w) - V/(z)) . (13) 
Taylor's formula up to order 2 for V/(w) around z yields 

l — l 2 

V/(w) = V/(z) + V 2 /(z)(w - z) + OC-—J-) , 

where V 2 /(z)(w - z) is the result of applying the matrix V 2 7(z) to the vector w - z. The 
mean value of V 2 I(z)(w - z) on B(z, t) is zero, by antisymmetry, and thus, by (fT3l) . 

|V ^ (Z)I - C 1^ 

and so, by (PT21) 

|£ f (x-y)-^(x)|<C|y|— . 
Integrating in ? we finally get 

/ f W/3 7 dtV \y\ I f w/3 \' M 



Interval 2: |x|/3 < t < 2\x\. The function I * Xt is continuously differentiable on 
R n \ 5 f , S t = {x : \x\ = t], because its distributional gradient is given by I * Vxt an d each 
component of Vxt is a Radon measure supported on 5,. The gradient of I * Xt is given at 
each point x e W \ S t by the principal value integral 



p.v.(VI *Xt){x) = -(n - \)p.v. j- — +1 



which exists for all such x. The difficulty in the interval under consideration is that it may 
happen that \x\ = t and then the gradient of I * Xt has a singularity at such an x. We need 
the following estimate. 



Lemma 1. 



< C log 



\x\ + t 
\x\-t\'- 



x e 



Proof. Assume without loss of generality that x = (jci, 0, . . . , 0). The coordinates y 7 , j ± 
1, change sign under reflection around the yi axes. Hence 



p.v 



J 

Jb(xJ) 



\y\ n+l 



dy = 0, 1 < j < n 
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Now, if \x I < t, 



If |x | > t, 



JB(x,t) 



\y\ 



n+ 1 



dy 



,,j 



B(x,t)\B(Q,Hx\) 



vln+l 



<c - 



Clog 



t + |jc| 
f - \x\ 



p- v - f n^r^ = f 

jb(xj) \y\ jb{ 



set,?) M 



;z+l 



dy 



< C 



r m+t dt _ 

J\x\-t t 



Clog 



|jc| + t 
x\ - 1 



Assume without loss of generality that y = (yi, 0, . . . , 0). The distributional gradient 
ofI*Xt is 

which is in L 2 . Then / *Xt e W 1 ' 2 and consequently is absolutely continuous on almost all 
lines parallel to the first axes. Therefore 



K t (x-y)-K t (x) = - VK t (x-ry)-ydT 



for almost all x and 



\y\ r l 



\K t (x-y)-K t (x)\<C^- 1+log 



x" 



\x - ry\ + t 
\x-ry\- t\ 



dr . 



Hence 

'2W 



/ p2\x\ 
VM/3 



\K t (x-y)-K t (x)\ 2 ^\ 



vl'l+l 



f p7M / s>l 
i Jw/3 \Jo 



1+log 



|x-ry| + f 
H* - ry| - ?| 



jc n+1 



where the last identity is a definition of D. Applying Schwarz to the inner integral in D 
and then changing the order of integration we get 



D 2 < 



M < nl\x\ 
10 WW/3 



1+log 



\x - Ty\ + t \ 2 dt^ 



\x-Ty\-t\) t ) 



dr. 
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For each r make the change of variables 



\x - ry\ 

to conclude that 



D 2 < | 1 1 + log 

2/9 



1 + s \ ds 



W-s\ 



Interval 3: 2\x \ < t. For each z in the segment [x - y,y] we have 5(0, t/4) c B(z, t). 
Then, by (HU), 



— r 



\B(z,t)\J BizJ) \ Bi0 , t /4 } \M n+l \z\ n+1 



and so 



|V^(z)|<C^-, ze[x-y,y] 
\x\" 



Hence, owing to (PT2l) . 



be 



and thus 



r°° 9 Lvi / r°° dtv \y\ 



>2\x\ ' / l A l \J2\x\ 

which completes the proof of the strengthened form of Hormander's condition (TTOl) . 

We turn now to prove that condition (2) in Theorem Q] is sufficient for / e W 1,p . Let 
/ 6 LP satisfy S\(f) e Z/. Take an infinitely differentiable function (f> > with compact 
support in 5(0, 1), f (f> = 1 and set <p £ (x) = pr^(f), e > 0. Consider the regularized 
functions f £ = f* <p £ . Then f £ is infinitely differentiable and ||V/ e || p < ||/|| p ||V0 e ||i, so that 
f £ £ W x,p . Thus, as we have shown before, 

l|V/ e ||-||Si(/ e )|l p . 

We want now to estimate \\S i(/ e )|| p independently of e. Since 

C/eW<) - M*) = df *Xt~f)* <Pe) (X) , 

Minkowsky's integral inequality gives 

Sl(f £ )(x) = \\(feh( x ,t) ~ /eWllLW) ^ ( S * , 

and so ||V/ e || < C ||5 i(/)|| p , e > 0. For an appropriate sequence 6j — > the sequences d k f £ . 
tend in the weak ★ topology of LP to some function gk e LP, 1 < < n. On the other 
hand, f £ — > f in LP as £ — > and thus d^/e — » dtf, 1 < & < n in the weak topology of 
distributions. Therefore d^f = gk for all k and so / e W x,p . 
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3 Proof of Theorem El 

The difficult direction is (1) implies (2) and this is the first we tackle. We start by showing 
that if / e W 2 ' 2 then 

IIW)ll2 = c||A/|| 2 (15) 

where the square function S 2 (f) is defined in ©. To apply Plancherel in the x variable it 
is convenient to write the innermost integrand in © as 



f ( 
"f ( 

JB(0,t) \ 



f(x + h)-f(x) 



(f 

\JB(x, t) 



In 



dz \\y - x\ dy 



j 

J 5(0,0 



Aftx + Q 
2n 



dk\\h\ 2 \dh. 



Applying Plancherel we get, for some dimensional constant c, 



c\\S 2 (f)\\ 2 2 = 



IJ 

Jr" Jb(o,o 



exp (igh) - 1 + 



(f 1 

\ J 5(0,0 



exp (igk) dk 



In 



/v 1 (XX 

dh\m\ 2 dt-. 



Make appropriate dilations in the integrals with respect to the variables h and k to bring 
the integrals on 5(0, 1). Then use that the Fourier transform of ^ 5(0,1) is a radial 
function, and thus of the form F(|£|) for a certain function F defined on [0, 00). The result 
is 



c\\S 2 (f)\\ 2 2 = f f 
Jr" Jo 

The change of variables r = 1 |£| yields 



Fim-i+^Fit^^-f \h\ 2 dh 



/B(0,1) 



CI I a o 

7 \m\ 2 d{. 



c||S 2 (/)||l = /||A/||2 



where / is the integral 



I 

Jo 



F(t) - 1 + t 2 F(t) 



Ln J 5(0,1) 



\h\ 2 dh 



dr 



The only task left is to prove that the above integral is finite. Now, as r — » 0, 

F( T ) - T exp(ihiT)dh 

J 6(0,1) 



"f ( 

J 5(0,1) \ 



1 + J^iT - ^l^" 2 + 



= 1 



z \ J 5(0,1 



h 2 dh\T 2 + 0(r 4 ). 



(16) 



11 



Hence 



F(r) - 1 + t 2 F(t)^- £ \h\ 2 dh 

2^ JB(0,1) 



\-\£ h 2 dh + ^-£ \h\ 2 dh) t 2 + 0(t 4 ) = 0(t 4 ) , 
\ 2 Jb(o,i) 2n Jb(o,i) / 



because clearly _£ (Q \h\ 2 dh = n f h 2 dh. Therefore the integral (PT6l) is convergent 
at t = 0. 

To deal with the case t -> oo we recall that F can be expressed in terms of Bessel 
functions. Concretely, one has ( llGrl Appendix B.5, p. 429]) 

P,1P(T) = ^. 

The asymptotic behaviour of Jn/iij) gives the inequality 

\F{r)\<C^<C 1 -, 



hit) = 77a ■ 



which shows that the integral (1161) is convergent at oo. 

We turn our attention to the case 1 < p < oo. Let I 2 (x) stand for the kernel defined on 
the Fourier transform side by 

1 

In other words, I 2 is minus the standard fundamental solution of the Laplacean. Thus 
I 2 (x) = c„ l/\x\ n ~ 2 if n > 3, 1 2 (x) = log \x\ if n = 2 and I 2 {x) = -\ \x\ if n = 1. Given 
any / 6 W 2 ' p there exists g e LP such that / = I 2 * g (indeed, g = -Af). We claim that 
there exists a singular integral operator T(g) taking values in L 2 (dt/t 5 ) such that 

S 2 (f)(x) = \\T(g)(x)\\ L 2 (dm . (17) 

Set ^ 

xix) = \mT)\ Xmj)(x) 

and 

1 x 
Xt(x) = -Xi~) ■ 



Then, letting M = f \z\ 2 dz, 



f ( 



B(0, 1) 

f(y) - fix) - ^(A/) e( ,, \y - *| 2 ) dy = ((/ 2 * Xt -I 2 -¥L t 2 Xt ) * g)(x) 

= (K t *g)(x), 
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where 

K t {x) = (I 2 * X ,){x) - I 2 (x) - t 2 Xt (x) . 

Setting T(g)(x) = (K t * g)(x) we get $TJ} from the definition of S 2 (f) in ©. Then O 
translates into 

r iii , gwig w /x=ciigii2, 

and we conclude that T is an operator mapping isometrically L 2 (R") into L 2 (R", L 2 (dt/t 5 )), 
modulo the constant c . If the kernel K t (x) of T satisfies Hormander's condition 

f - J) - K t (x)\\ LHdt/tS) <C, ye R", 

J\x\>2\y\ 

then a well known result of Benedek, Calderon and Panzone on vector valued singular 
integrals (see HGRl Theorem 3.4, p. 492]) yields the LP estimate 

which can be rewritten as 

\\S2{f)\\ P <c P w\\ P - 

The reverse inequality follows from polarization from (TT5T) by a well known duality argu- 
ment ( HGR1 p. 507]) and so the proof of the necessary condition is complete. 

We are going to prove the following stronger version of Hormander's condition 

\\K t (x - y) - K t (x)\\ L 2 {dtlfi) < C \x\ > 2\y\ . (18) 

\x\ 

For this we deal separately with the kernels H,(x) = (I 2 *Xt)( x ) ~ h(x) and t 2 Xt(x). For 
t 2 Xt{x) we first remark that the quantity \xt(x-y) ~Xt(x)\ is non-zero only if \x-y\ < t <\x\ 
or \x\ < t < |x-y|, in which cases takes the value 1 /c„ t n , c„ = \B(0, 1)|. On the other hand, 
if |x| > 2\y\ then each z in the segment joining x and x — y satisfies |z| > \x\/2. Assume that 
|x - y\ < \x\ (the case |x| < \x — y\ is similar). Then 




We check now that H t satisfies the stronger form of Hormander's condition. If t < 
\x |/2, then the origin does not belong to the ball B(x - y, t) nor to the ball B(x, f). Since 
I 2 is harmonic off the origin, the mean of I 2 on these balls is the value of I 2 at the center. 
Therefore H t (x - y) - H t (x) = in this case. 



13 



If t > |x|/2, then 

\H t {x-y)-H t {x)\<\y\ sup \VH t (z)\ < C . 
The last inequality follows from 

V# f (z) = £ V/ 2 (w) Jw - V/ 2 (z) , 

|V/ 2 (z)| < C 1/lzr 1 < C and 

| £ VI 2 (w)dw\ < £ — !— < C -I— . 

Therefore 

Jo ur 1 Uw/2 w n+1 

We turn now to prove that condition (2) in Theorem [2] is sufficient for / e W 2 ' p . 
Let / and g in V satisfy S 2 (/, g) s Z/. Take an infinitely differentiable function (p > 
with compact support in 5(0, 1), f (/> = 1 and set (f> € (x) = pr^(f), e > 0. Consider the 
regularized functions f e = f * (p e and g e = g * e . Then f e is infinitely differentiable and 
l|A/ e ||p < ll/ll P l|A0 e ||i, so that f e 6 IV 2 -". Recalling that M = f B{0l) \z\ 2 dz, we get, by 
Minkowsky's integral inequality, 

S 2 (fe,ge)(x) = \\(f £ * X tXx) ~ MX) ~ (g e *Xt){x) M 2 t 2 \\ L 2 (dt/t5) 

= II (if *Xt)-f~(g *Xt) M 2 1 2 * <f> e ) W|| L 2 W , 5) 
<(S 2 (f,g)*cf> e )(x). 

Now we want to compare (1 /2n)A/ e and g £ . Define 

2 1/2 



M 2 



^-(A/ e *^)(jc) - (g 6 *xt){x) 
2n 



dt 



Then 

D e (x)<S a (/ e )(x) + S a (/ e ,g e )(x) 

<S a (f £ )(x) + (S a (f,g)(x)*fc)(x), 
and thus Z) e (x) is an LP function. In particular D £ (x) < oo, for almost all jteR". Hence 

|(l/2n)A/ e (x) -g e (x)\ = lim|(l/2n)(A/ e **,)(*) - (g e **,)(*)! =0, 

/->o 

for almost all x e I", and so (l/2«)A/ e — > g in LP as e — > 0. Since / e — > / in Z/ as 6 — » 0, 
then A/ e — » A/ in the weak topology of distributions. Therefore (l/2n)A/ = g and the 
proof is complete. 
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4 Proof of Theorem ffl 



The difficult direction in Theorem[3]is to show that condition (2) is necessary for / 6 W a ' p . 
The proof follows the pattern already described in the preceding sections. One introduces 
an operator T taking values in L 2 (dt/t 2a+l ) and shows via a Fourier transform estimate that 
T sends L 2 (R") into L 2 (R.",L 2 (dt/t 2a+l )) isometrically (modulo a multiplicative constant). 
The second step consists in showing that its kernel satisfies Hormander's condition, after 
which one appeals to a well known result of Benedek, Calderon and Panzone on vector 
valued singular integrals to finish the proof. 



4.1 The fundamental solution of (-A) a/2 

Let I a be the fundamental solution of (-A) ff/2 , that is, I a is a function such that I a {^) = \£\~ a 
and is normalized prescribing some behavior at oo. It is crucial for our proof to have an 
explicit expression for I a . The result is as follows (see HACL1 or HMOPVl p. 3699]). 
If a is not integer then 

/«(*) = c a , n \x\ a - n , xeR n , (19) 

for some constant c a>n depending only on a and n. 

The same formula works if a is an even integer and the dimension is odd or if a is an 
odd integer and the dimension is even. 

The remaining cases, that is, a and n are even integers or a and n are odd integers are 
special in some cases. If a < n formula (fT9l still holds, but if a is of the form n + IN, for 
some non-negative integer N, then 

I a (x) = c a>n \x\ a ~ n (A + B log \x |), x e R" , 

where c aA , A and B are constants depending on a and n, and B £ 0. Thus in this cases a 
logarithmic factor is present. For instance, if a = n, then I a (x) = B\og\x\. If n = 1 and 
a = 2, then I 2 (x) = -(1/2) \x | and there is no logarithmic factor. 



4.2 The case p = 2 

Given a positive real number a let be the unique integer satisfying 2N < a < IN + 2. 
Define the square function associated with a by 



Sa(f)\x) 



= r f 

Jo JB(x,t) 



PN(y, x) 



dy 



dt 

— , X G 



(20) 



where p N (y, x) is 



/Cy)-/(*)-^A/(;t) \y-x\ 2 - 
2n 



^-A^-V(x) \y - x\ 2 ^ - ^-(A N f) Bix Jy - x\ 2N , 
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Recall that Lj = A ; (|x \ 2j ) and that the role which the Lj play in Taylor expansions was 
discussed just before the statement of Theorem [3] in the introduction. 
In this subsection we prove that 

l|5 ff (/)||2 = c||(-Ar /2 (/)|| 2 . (21) 

Our plan is to integrate in x in (1201) . interchange the integration in x and t and then 
apply Plancherel in x. Before we remark that making the change of variables y = x + th 
we transform integrals on B(x, t) in integrals on B(0, 1) and we get 

f p N (y, x)dy= f f(x + th) dh - V t 2J f \h\ 2j dh 

JB(x,t) Jfi(O.i) L i Jb(0,\) 

- f A N f(x + th) dh t 2N f \h\ 2N dh . 

Jb(0, 1) J 5(0,1) 

Now apply Plancherel in x, as explained before, and make the change of variables r = t \%\, 
where £ is the variable in the frequency side. We obtain 

\\S a (f)\\ 2 2 = cI\\(-A) a/2 f\\ 2 2 , 

where 

/ = f F(r) - Y(-l) j T 2 ^ - (-l)V^ F(r) 
Jo ^ L j l n 

Here F is a function defined on [0, <x>) such that F(\g\) gives the Fourier transform of 
the radial function ^ 8(0,1) a * the point and we have introduced the notation My = 

f j |/z| 2i d/z. We have to show that the integral / is finite. 

Using the series expansion of the exponential we see that, as r — > 0, 

F(t) = -r exp («/zit) rf/i 

Jfi(0,l) 



2 

r 2a+l 



= 1 + ■ • ■ + (-If r 2N -L- f *f <ft + 0(t 2N + 2 ) . 

(Z-/Vj! Jfi(0,l) 

We need to compare f B(Q hf dh with f m 1} |/*| 2W d/z. The linear functionals P -> A 2 ;'(P) 

and P — > J^ (Q P, defined on the space // 2 j of homogeneous polynomials of degree 2 _/', 
have the same kernel. This follows from the discussion before the statement of Theorem^] 
in the introduction. Therefore, for some constant c, 



A 2 ^(P) = c f 



P P e H 2j 

B(0,1) 



Taking P(x) = \x \ 2j we get Lj = c f 1} \ x \ 2j = dx, and taking P{x) = x 2j we get 
(2 ^ ! = c £(o,i) x i Jrfx Hence 



— ^— x 2j dx= -j- -f \x \ 2j dx = — - 

(2j)- Je(o,i) Lj Js(o,i) Lj 
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and thus, owing to the definition of I and the fact that F(r) = 1 + 0(t 2 ), as r — > 0, 

dj 



o 



6>(r 2(2/v+2) ) 



r 2«+l 



as r — » . 



Then I is convergent at because a < 2N + 2. 

We turn now to the case r — > oo. Notice that the only difficulty is the last term in the 
integrand of I, because 



x>0 

Ji 



A' 



dr 



'— T <oo, 0<j<iV-l, 
provided 2/V < a. To deal with the term 



Ji 



r 2w F(r)r 



r 2a+l 



(22) 



we only need to recall F can be expressed in terms of Bessel functions. Concretely, one 
has (IS Appendix B.5, p. 429]) 



|5(0, 1)| F(t) = 



•In/2 



The asymptotic behaviour of /„/2(t) gives the inequality, as r — > oo, 

|F(r)|<c4r^, 
which shows that the integral (1221) is convergent finite provided 2/V < a. 



4.3 A vector valued operator and its kernel 

Given / e W ^, there exists a function g e Z/ such that / = I a *g. Indeed, g = (-A) a/2 (f). 
Then 



r ( N-l 1 

f /(y) - /(*) - X y- A j f(x) \y - x\ 2j - -(A7) B( , () \ y - x\ 2N 

JB(x,t){ j^O L J LN ) 



dy = (K t *g)(x), 



where the kernel K t (x) is 

/ N-l 



K t (x) = f I a (y) -Yj- &I a (x) \y - x\ 2j - ^-(A N I a ) B(x , t) \y - x\ 2N 

JB(x,t)\ j^O J N 

Hence the square function associated with the smoothness index a is 

dt 



dy . (23) 



S a (f) 2 (x)= \(K t *g)(x)\ 2 
Jo 



X € 
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Define an operator T acting on functions / e L 2 (R") by 

Tg(x) = (K t *g)(x), xeR". 

The identity dZD in subsection 4.2 says that T takes values in L 2 (R n ,L 2 (dt/t 2a+l )) and, 
more precisely, that 



IIWfe(A/^)^=ll' S «COIll = c||g|ll 



Therefore T is an operator mapping isometrically (modulo a multiplicative constant) 
L 2 (R") into L 2 (R", L 2 (dt/t 2a+i )) and we have an explicit expression for its kernel. If we 
can prove that K t (x) satisfies Hormander's condition 

f ||* f 0c - y) - K t (x)\\ L 2 (dt/t 2 a+i) dx<C, y e R" , 

then the proof is finished by appealing to a well known result of Benedek, Calderon and 
Panzone ( HGRl Theorem 3.4, p. 492]; see also HGR1 p. 507]). In fact, we will show the 
following stronger version of Hormander's condition 

\\K t (x - y) - K t (x)\\ LHdt/t 2 a+l) < C \x\ > 2\y\ , (24) 

l-^l 

for some y > depending on a and n. 

The proof of (12"41 is lengthy. In the next subsection we will consider the case of small 
"increments" in t, namely t < \x\/3. 

4.4 Hormander's condition: t < \x\/2> 

We distinguish two cases: IN < a < IN + 2 and a = IN. To deal with the first case we 
set K t (x) = K\ l \x) - Kf\x), where 

N , \ 

dy (25) 



K\'\x) = f I a (y) - Y 1 A j I a (x) \y - x\ 2j 

Jb(xJ){ 



and 

Kf\x) = ( -3-( (A N I a ) B(x , t) - A N I a (x) ) \y - x\ 2N dy . 

JB{x,t) L N 

We first estimate To compute the gradient of K$ we remark that 

N 



K\ l \x) = (4 * Xt){x) - Y ^ A'/ a (x) , 



,/ : 



=0 
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where Mj = f \z\ 2j dz. Thus 



JB(0,l) 

r ( N 1 "1 

VK^ix) = f VI a (y) - V - A^'(V4)(x) |y - x| 2 ' Jy. 

Let P ra (F, x) stand for the Taylor polynomial of degree m of the function F around the 
point x. Therefore 



VK { t l \x) = f <yi a (y) - P 2N+ i(VI a , x)(y)) dy , 

)B(x, t) 

because the terms which have been added have zero integral on the ball B(x, t), either 
because they are Taylor homogeneous polynomials of VI a of odd degree or because they 
are the "zero integral part" of a Taylor homogeneous polynomial of VI a of even degree 
(see the discussion before the statement of Theorem 3 in the introduction). Given x and y 
such that |x| > 2\y\, apply the formula above to estimate V^ (1) (z) for z in the segment 
from x - y toy. The standard estimate for the Taylor remainder gives 

\VK ( t l \z)\<t 2N+2 sup |V 2/v+3 4(w)| . 

weB(z, f) 

Notice that if z e [x - y,y], w e B(z, t) and t < |x|/3, then \w\ > |x|/6. Now, one has to 
observe that 

\v 2N+3 i a (w)\ < c \ w r n ~ 2N - 3 , 

owing to the fact that possible logarithmic terms do not appear because the exponent a - 
n - IN - 3 < —n — 1 is negative. By the mean value Theorem we then get 

Kl)/^ „\ vW/^W ^ U,l \T7rWf„\\ ^ r U,l *2W+2 i \a-n-2N-3 



^B( 



\K) l \x - y) - K^(x)\ < \y\ sup |V£^(z)| < C |y| t lN+l \x\ 



Since 



x|/3 dt \ m 



we obtain 

w/3 dt \ 1/2 Ivl 

n I^U-^'W^) SCjM.. (26) 

which is the stronger form of Hormander's condition (1241) with 7= 1 in the domain t< |x|/3. 
We turn now to estimate K (2) . Arguing as above we get 



VK^ 2 \x) = C f (A N (VI a )(y) - J P 1 (A /V (V4), x)(yj) \y - x\ 2N dy . 

JB(x, t) 

If z e [x-y,y], then 

\VK?\z)\ < t 2N + 2 sup IV^/^Cw)! < Cf 27V+2 ^1—2^-3 ^ 

weB(z,t) 

and so we get (T26l) with K ( t l) replaced by K (2) exactly as before. 
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Let us consider now the case a = 2N. Since A n I 2 n is the Dirac delta at 0, A N I 2N (x) = 0. 
Hence K t (x) = K\ l \x) - K\'(x), where is given by (1251) with a replaced by 2N and 

K?\x) =f f( (A N I 2N ) B(xJ) \y - x\ 2N dy = ^ t 2N (A N I 2N ) B(xJ) . 

The kernel K\ l) is estimated as in the first case by just setting a = 2N. The kernel K t 
requires a different argument. 
Set 

1 

x{x) ~ \mr)\ XBm){x) 

and 

1 x 
Xt(x) = -X(-) , 

so that, since A N I 2N is the Dirac delta at the origin, K (2) is a constant multiple of t 2N Xt- 
We show now that this kernel satisfies the strong form of Hormander's condition. The 
quantity \x t (x - y) - Xt(*)\ 1S non-zero only if \x - y\ < t < \x\ or \x\ < t < \x - y\, in 
which cases takes the value l/c„ t'\ c„ = \B(0, 1)|. On the other hand, if \x\ > 2\y\ then 
each z in the segment joining x and x-y satisfies \z\ > \x\/2. Assume that \x-y\ < \x\ (the 
case |x| < \x — y\ is similar). Then 



< c 13 



\x - y\ 2n \x\ 2n ) \x\ n+l12 ' 
which is (|2D with y = 1/2. 

4.5 Hormander's condition: t > \x\/3 

The last term in the definition (1231) of the kernel K, is of the form 

- f ±(A N I a ) B(x , t) \y - x\ 2N dy = t 2N (A N I a ) B(x , t) , 

JB(x,t) 

and our first goal is to show that this kernel satisfies Hormander's condition in the do- 
main t > \x |/3. Notice that 

provided no logarithmic factors appear, which is the case if a - 2N < n. Since a - 2N < 2 
the inequality a - 2N < n certainly holds whenever n > 2. Hence the only cases with 
logarithmic factors are n = 1 and a = 2N + 1, which will be dealt with separately. 
In the previous subsection we dealt with the case a = 2N and so we can assume that 
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IN < a < 2N + 2 , which implies that A N I a is locally integrable. Given x and y with 
I* I > 2\y\ let D stand for the symmetric difference (B(x-y, t)\B(x, t))V(B(x, t)\B(x-y, t)). 
It is easy to realize that \D\ < C |y| f 1 ' 1 . We remind the reader of the following well known 
and easy to prove inequality. 

Lemma. Let E be a measurable subset ofW and < J3 < n. Then 

dz 



< C \Ef /n , 



where \E\ is the Lebesgue measure ofE. 
We have 



" f- 

J D \Z\"- 



r\(A"I a ) Bix - y , t) -(A"I a ) B(x , t) \<Cf 

< gj2N-n^f-iy.a-ZN)/n _ ^(a-ZN)/n t <i-n-(a-2N)/n 



if 

\J\x\i 



and 

dt \ m 

f 2(a-n-(a-2N)/n) ul \ „ i i-n-(a-2N)/n 
'W/3 1 

Combining the preceding two inequalities one gets Hormander's condition (1241) with y = 
(a - 2N)/n in the domain t >\x |/3. 

We have to investigate now the exceptional cases n = 1 and a = 2N + 1 , in which no 
logarithmic factors appear. We have 

I a (x) = \x\ 2N (A + B\og\x\), 

for some constants A and 5^0. Then 

A N I a (x) = (A' + B'\og\x\) 

and 

-fA N I a (x) = Cv.p.-, 
ax x 

in the distributions sense. Thus, denoting by H the Hilbert transform and applying 
Schwarz inequality, 

t 2N \(A N I a ) B{x ^ t) - (A N I a ) B(x J = t 2N \ f H(xt)(j) dr\ 

J x-y 

< t 2N \y\ l/2 ||*|b = t 2 ^' 2 \y\" 2 . 

Since 



If 

\J\x 



t 2(2N-m)jL) m = C \x\ 2N ~ ll2 - a = C \xV l ~ m 



,2a+l 

x\/3 I 
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we obtain Hormander's condition d2U) with y = 1/2 in the domain t > \x |/3. 

Our plan is to continue the proof by distinguishing three cases: a < n + I, a = n + 1 
and a > n + 1 . 

If a < n + 1, then all terms in the expression (1231) defining K, satisfy Hormander's 
condition in the domain t > \x |/3. Indeed, consider first the terms of the form t 2j A->I a (x), 
j > 0. One has the gradient estimate 



\t 2j VA j I a (x)\ < C t 2j \xr n ' 2j ' 1 , 



(27) 



because no logarithmic factors appear, the reason being that the exponent a -n-2j - \ < 
a - (n + 1) is negative. Since 



( r ^•)^ r ) 1/2 

VM/3 t 2a+l ) 



= C \x\ 2j ~ a , 



we get Hormander's condition with y = 1 in the domain t >\x |/3. 

It remains to look at the first term JF f) I a (y)dy in (1231) . Its gradient can be easily 
estimated as follows 



If ' 

\JB(x,t) 



< C 



f \yr n ~ 

JB(x,t) 



dy . 



Notice that there are no logarithmic factors precisely because a < n + 1 . The integrand 
in the last integral is locally integrable if and only if a > 1. Assume for the moment that 
1 < a < n + 1 . Then 

f VI a (y)dy < Cr n+a - 1 . 

>B(x, t) 



Since 



if 

WM/3 



t 2(a-n-\) 



dt 

t2a+l 



1/2 



C \x 



-n-1 



we get Hormander's condition with y = 1 in the domain t > \x |/3. The case a = 1 has 
been treated in section 1, so we can assume that < a < 1. In this case, denoting by D 
the symmetric difference between B(x, t) and B(x - y, t), we obtain 



j 

Jb< 



I a (y) dy-f I a (y) dy\ < C t~" f 

B(x-y,t) jb(xj) Jd \y\ 



Since 



(f 

VM/3 



f 2(cY-n—a/ri) 



<Ct" (t n ~ l \y\) a/n = Ct a ~"~ a/n \y\ a/n 
dt 



1/2 



= C\x 



—n—a/n 



we get Hormander's condition with y = a/n in the domain t >\x |/3. 

We tackle now the case a = n + 1. Since a and n are integers with different parity no 
logarithmic factor will appear in I a . Thus I a (x) = C\x\. The proof above shows that the 
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terms t 2 -> A j I a (x) appearing in the expression (1231) of the kernel K t still work for j > 1 . 
The remaining term is 

(I a (y)-I a (x))dy 



^Bt 



B(x, t) 

and its gradient is estimated by remarking that the function \x | satisfies a Lipschitz condi- 
tion. We obtain 



But clearly 



f 1 

JB(x, t) 



lXl/3 *H 



(V4(y) - VI a (x)) dy 



<C 



1/2 



C \x\~ a = C \x\ 



-n-l 



which completes the argument. 

We turn our attention to the case a > n + 1 . There is a unique positive integer M such 
that -1 < a — n — 2M < 1 . The part of K, which has to be estimated is 



H t (x) = f I a (y) -V- A j I a (x) \y-x 

JB(x,t){ 



2j 



dy . 



We split H t into two terms according to M, that is, H t = // f (1) - Hf\ where 

M-l 



// f (1) w 



= f /«(y)- X -5-A'4(x)|y-x| 



2,/ 



and 



Hf\x)= \ Y-^I a {x)\y-x\ 2j 

JB( X J){ j=M Lj j 



V/ 
,(2) 



dy . 



The estimate of each of the terms in H) is performed as we did for the case a < n + 1 . 
The gradient estimate is exactly (1271) . Now no logarithmic factors appear because the 
exponent a -n-2j -1 < a - n - 2M - 1 < is negative. The rest is as before. 

To estimate H { t Y) we distinguish three cases: -1 < a - n - 2M <0, < a-n- 2M < 1 



and a-n - 2M = 0. In the first case we write the gradient of H) as 

r ( M ~ l i 

VH^ix) = \ VI a (y) -V- A'(V/ a )(x) |y - x\ 2j dy 

Jb(xj){ ) 



j 

JB(xJ) 



(V4(y)-P 2 M- 2 (V4,x)(y)) dy, 



where P2M-2 is the Taylor polynomial of degree 2M - 2 of V/ a around the point x. As 
before, the added terms have zero integral on B(x, t) either because they are homogeneous 
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Taylor polynomials of odd degree or the "zero integral part" of homogeneous Taylor poly- 
nomials of even degree. Now fix y in B(x, t) but not in the half line issuing from x and 
passing through the origin. Define a function g on the interval [0, 1] by 

g(r) = VI a (x + r(y- x)) - P m . z <yi a , x)(x + r(y- x)), < r < 1 . 

Since g j) (0) = 0, < j < 2M - 2, 



2M-2 



V/„(y) - PiM-iWIa, x)(y) = g(l) - ^ 



g j) (0) 



(l-T) 



2M-2 



g 2M - l \r)dr, 



o (2M-2)! 



by the integral form of Taylor's remainder. The obvious estimate for the derivative of g of 
order 2M - 1 is 



|g 2M -D( T )| < IV 2 ^ 1 V/ tt (jc + r(y - x))\\y - x\ 2M - y < C 



t 2M-l 



\x + r(y - x )\"^ a ' 2M ^ 

Since we are in the first case a is not integer and thus no logarithmic factors exist. More- 
over < n - (a - 2M) < 1, which implies that and that 1 /\z |' i_(Q,_2M) is locally integrable 
in any dimension. Therefore 



= Ct 



|V# f (1) (jc)| < C f (t 2M - l - n f ■ 

Jo \ JB(x,t) \ x + v 

J0 \JB(x,tT 

Jo 

_ f-n-\ C ^_ 

JO T 



dy 



(y - jc)|"-( ff - 2M > 



dr 



n-(a~2M) I T n 



2M-\-n 



{tTT -2Mdr 



Since 



if 

WM/3 



f 2{a-n-Y) 



dt 

t2a+\ 



1/2 



C\x 



-n-l 



we get Hormander's condition with y = 1 in the domain t >\x |/3. 

Let us consider the second case: < a-n-2M < 1 . This time we express the gradient 
of// f (1) by means of a Taylor polynomial of degree 2M - 1 : 

r ( M ~ l i 

VH^ix) = f V4(y) - V - A\VI a )(x) \y - x\ 2j dy 
Jb(xj){ J 



JB(x,t) 



(V4Cy)-P 2M _ 1 (V/ ff ,x)(y)) dy. 
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Using again the integral form of the Taylor remainder of the function g, with P2M-2 re- 
placed by P2M-1, we obtain 

"Ilf 

J() \JB(xJt) 



ct 2M - 



dz \ dr 

^ |^ |n-(a-2M-l) j T n 



<C7 2M "« f 1 (/T)^ 24f - 1 — 

Jo 

= ^-»- 1 f 1 — = cr 

J T n-(or-2M-l) 



.or— n— 1 



from which we get the desired estimate as before. 

We turn now to the last case left, a = n + 2M, with M a positive integer. In this case 

I a (x) = C\x\ 2M (A + B\og\x\), jceR", fl^O, 

where A, B and C are constants depending on n and M. We also have 

A m_1 4(jc) = C \x\ 2 (A 1 + B x log \x |), jc e W 

and 

VA M ~ 1 4(x) = Cx(A 2 + 5 2 log \x |), xel". 
In particular VA^ -1 /^ is in the Zygmund class on K". We have 

M-2 . . \ 



r 1 1 1 

V// r (1) (x) = f V/ a (y) - V - A^'(V4)(x) ly-xp - A^V/^*) ly-*! 2 ^ 2 

= f |V/ (y) - - 

^B(.v,f)\ 



P 2M -3(V/ a , x)(y) - A^W/JOO |y - xl^ 2 | . 

Introduce the function g as above, with P2M-2 replaced by P 2 m~3, so that 



2M-3 



V4(y) - P 2 m-3(V4, *)(y) = g(l) - 2 



'(0) 



j=o J 



= (2M-2)(l-r) 2M - 3 5,, \;, Jr. 
Jo 



(2M-2)! 
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Now 



g 2M - 2 \r) 
(2M-2)! 



E 



=2M-2 



dPVI a (x + r(y - s)) 



(y-xf 



(y-xf 



= Yj (^V4(x + r(y-x))-^V4(x)) 

l6|=2M-2 



=2M-2 



The last term in the preceding equation is the homogeneous Taylor polynomial of de- 
gree 2M - 2 of the vector VI a around the point x. It is then equal to a homogeneous poly- 
nomial of the same degree with zero integral on B(x, t) plus A M ~ l (VI a )(x) \y - x\ 2M ~ 2 
(by the discussion before the statement of Theorem 3 in the introduction). Hence 



JB(x,t) 





(dPVI a (x)\ 


E 

\]p\=aM 


-2' ^ ! / 



(y-xf 



1 



A M -\VI a )(x) \y-x 



2M-2 



dy = 0, 



and therefore, remarking that £(2M - 2)(1 - T ) 2M ~ 3 <ir = 1, 



VH^ix) 



Jb(xj) Jo 



(2M-2)(1-t) 



2M-3 



£ (^V4(jc+r(y-x)) - ^V/ a (x)) 



VP=2M-2 



Thus 



|V// f (1) (x)| < C f V f (^VUx + t^-i))-^VUx)) 



(y-*y 



dr dy . 



dy 



dr . 



Making the change of variables h = r(y - x) the integral in y above becomes 

f hP 

j = T -w \ (dPm a (x + h)-dPvi a (x))-dh, 

JB(0,tr) P- 

which is invariant under the change of variables h' = -h, because \B\ is even. Hence 

hP 



2J = T 



\ (dPVI a (x + h) + dPVI a (x -h)-2 dPVI a (x)) — dh . 

Jb(0jt) P- 



Now we claim that dPVI a is in the Zygmund class for \B\ = 2M - 2. This follows from the 
fact that the Zygmund class in invariant under homogeneous smooth Calderon -Zygmund 
operators, A M ~ l is an elliptic operator and A M ~ l VI a is in the Zygmund class. Hence 



\J\<Ct~ w f \h\ l+ ^dh<C 

JB(0jt) 
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Thus 

IVfl^Oc)! <Ct 2M ~ l . 

Since 

we get Hormander's condition with y = 1 in the domain t >\x |/3. 

4.6 The sufficient condition 

In this section we prove that condition (2) in Theorem [3] is sufficient for / 6 W a,p . Let 
f,gi, . . . ,gN e L p satisfy S a (f,gi, ■ ■ ■ ,gN) e L p . Take an infinitely differentiable func- 
tion (p > with compact support in 5(0, 1), f (p = 1 and set <p e (x) = pr^(f), e > 0. Con- 
sider the regularized functions f £ = f* 4> e , gj, £ = gj *(p e , I < j < N. We want to show first 
that the infinitely differentiable function f e is in W a ' p . We have (-A) a/2 f £ = f * (-A) a/2 (p £ . 
We need a lemma. 

Lemma 2. 

(i) If <p is a function in the Scharwtz class and a any positive number, then (-A) a ^ 2 (p 
belongs to all L q spaces , 1 < q < oo. 

(ii) Iff e LP , 1 < p < oo , then (-A) a/2 f is a tempered distribution. 

Proof. Set i// = (-A) a/2 <p. If a = 2m with m a positive integer, then i// = (-A) m (p is in the 
Scharwtz class and so the conclusion in (i) follows. If a = 2m + 1, then 

n 

if, = (-A) 1/2 (-A)"V = -iJ] Rj( dj (-A)'V ) , 

7=1 

where Rj are the Riesz transforms, that is, the Calderon-Zygmund operators whose 
Fourier multiplier is It is clear from the formula above that i// is infinitely differ- 

entiable on R" and so the integrability issue is only at oo. Since dj (-A) m (p ), has zero 
integral, one has, as x — > oo, < C and so the conclusion follows. 

Assume now that m - 1 < a < m, for some positive integer m. Thus 

If m is even, of the form m = 2M for some positive integer M, then 

ifr = A M <p * I m . a , 

where I m - a (x) = C \x\'"~ a ~". Hence iff is infinitely differentiable on R". Since A M (p has zero 
integral, |i^(jc)| < C \x\ n ~ a ~"~ l , as x — » oo. But a-m+l > and thus <// is in all L q spaces. 
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If m is odd, of the form m = 2M + 1 for some non-negative integer M, then 

n 

i{, = -iJ]Rj(djA M ( p)*I, n ^. 

7=1 

Again iff is infinitely differentiable on W and, since Rj(djA M (p) has zero integral (just 
look at the Fourier transform and remark that it vanishes at the origin), we get |t//(x)| < 
C \x\ m ~ a ~"~ l , as x — > oo, which completes the proof of (i). 

To prove (ii) take a function </? in the Schwartz class. Let q be the exponent conjugate 
to p. Define the action of (-A) a/2 f on the Schwartz function <p as (f, (-A) a/2 cp) . By part 
(i) and Holder's inequality one has 

\((-Af l2 f,ip)\ = \(f,(-A) a/2 <p)\ <C\\f\\ p \\(-Ar /2 <p\\ q , 

which completes the proof of (ii). □ 

Let us continue the proof of the sufficiency of condition (2). By the lemma (-A) al2 (j) £ 
is in L 1 and so 



-AT /2 f £ \\ P = \\f*(-A) a '^ £ \\ p < IK-Ar^iii 



\a/2 



a/2. 



Hence f £ 6 W a ' p . 

Next, we claim that 



Sa(f £ ,gl,e, ■ ..,g N , £ )(x) < (S a (f,g U . . . , g N ) * (p £ )(x), 



(28) 



One has 



N-l 



S a (f, gu . . . , g N )(x) = ||(/ **,)(*) - f(x) - ^ Mj gj (x) & - M N (g N **,)(*) t 2N \ 



\L 2 {dt/t 2a+i ) 



7=1 



where Mj = j \z\ 2] dz. Minkowsky's integral inequality now readily yields (1281) . 



Set 



N-l 



D £ (x) = \\Y J M j ( 

7=1 



- gj, £ (x)) t 2J - M N I (— - - g N , £ ) * X t I (x) OIlW" 1 ) 



L 



'/V 



By 

D £ (x) < S a (f £ )(x) + S a (f £ ,gl, e, ■ ■ .,g N ,e)(x) 

< S a (f £ )(x) + (S a (f, gi,...,g N )* <p £ ){x) , 
and so D £ e LP. In particular, D £ (x) is finite for almost all x e W. Thus 

N-l 



lim inf 

f->0 



g Mj (^^ - gjAx)) P - M N ((^ - g N , £ ) (x) t 2N 



r = o, 
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for almost all x e R". It is easy to conclude that the only way this may happen is whenever 

A j f e (x) 



g j>£ (x), l<j<N. 



for almost all x e R n . Hence 

A j f e 



gj, l<j<N, 



in LP as e — > 0. Since / e — > / in Z/ as 6 — » 0, 

A'/ e ->A'/, l<j<N, 
in the weak topology of tempered distributions. Hence 

tJf 



Lj 



gj, l<j<N. 



We claim now that the functions f e are uniformly bounded in W a,p . Indeed, by the 
proof of necessity of condition (2) and by (1281) . 

||(-Ar /2 / e |l P - \\S a (f £ ,AfJL u ...,A N f £ /L N )\\ p 

<\\S a (f,Af/L u ...,A N f/L N )\\ p <oo. 

Hence there exist a function h e LP and a sequence e, — » as j — * oo such that 

(-Af/ 2 f £ .^h as j^oo 

in the weak ★ topology of V . On the other hand, by Lemma 2, {-A) al2 f is a tempered 
distribution and so 

{-Af l2 f e -> (-A) ff/2 / as 6^0 

in the weak topology of tempered distributions. Therefore {-A) al2 f = h £ L p and the 
proof is complete. 



5 Final remarks 



Let (X, d,fi) be a metric measure space, that is, X is a metric space with distance d and yu 
is a Borel measure on X. We assume that the support of \x is X. Then, given a > and 
1 < p < oo, we can define the Sobolev space W a ' p (X) as follows. Let Af be the unique 
integer such that 2Af < a < IN + 2. Given locally integrable functions /, gi, . . . , g N define 
a square function by 



S a (f,gl,g2,---,gN)(xY 



-flf 

JO JB(x,t) 



R N (y, x) 



dy.(y) 



dt 
7' 



X € 
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where D is the diameter of X and Ru(y, x) is 
R N (y, x) = f(y) - f(x) - 8l (x) d(y, x) 2 + ---- g^x) d(y, x) 2(N ~ l) - (g N ) B{x , t) d(y, x) 2N . 

Here the barred integral stands for the mean with respect to p. on the indicated set, B(x, t) 
is the open ball with center x and radius t and g B ( X ,t) is the mean of the function g on B(x, t). 

We say that a function / belongs to the Sobolev space W a,p (X) provided / e L p (jj) 
and there exist functions g u g 2 , . . . ,g N € L p (p) such that S a (f, gugz, ...,g N )& L p (pi). 

We have seen in the previous sections that this definition yields the usual Sobolev 
spaces if X = R" is endowed with the Euclidean distance and p is Lebesgue measure. One 
can prove with some effort that the same is true if R" is replaced by a half-space. Very 
likely this should also work for smoothly bounded domains, but we have not gone that 
far. 

There are many interesting questions one may ask about these new Sobolev spaces. 
For instance, how do they compare, for a = 1, to the known first order Sobolev spaces, 
notably those introduced by Hajlasz in [|H] or the Newtonian spaces of [[Sj ? For higher 
orders of smoothness one would like to compare them with those introduced by Liu, Lu 
and Wheeden in ULLWL One may also wonder about their intrinsic properties, namely, 
about versions of the Sobolev imbedding theorem, the Poincare inequality and so on. 

For the Sobolev imbedding theorem the following remark might be useful. In R" 
the L p space can be characterized by means of the following "zero smoothness" square 
function: 

So(f)\x)= |/ B( ^-/ BOc>2 o| 2 -, xgR". 
Jo t 

The result is then that a locally integrable function / is in U if and only if S (f) e U . The 

proof follows the pattern described several times in this paper. One first deals with the 

case p = 2 via a Fourier transform computation. Then one introduces a L 2 (dt/t)-\alued 

operator T such that 

\\T(f)\\ L 2 m „ c\\So(f)h 
and one shows that its kernel satisfies Hormander's condition. 
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A new characterization of Sobolev spaces on W 

Roc Alabern, Joan Mateu and Joan Verdera 



o 

CN ! Abstract 

q ■ In this paper we present a new characterization of Sobolev spaces on MP. Our 

characterizing condition is obtained via a quadratic multiscale expression which ex- 
ploits the particular symmetry properties of Euclidean space. An interesting feature 
■ of our condition is that depends only on the metric of R" and the Lebesgue measure, 

so that one can define Sobolev spaces of any order of smoothness on any metric 
^ , measure space. 

u 

"5 : 1 Introduction 

In this paper we present a new characterization of the Sobolev spaces W a ' p on R", where 
the smoothness index a is any positive real number and 1 < p < oo. Thus W a ' p consists 
of those functions / 6 LP = L P (R") such that (-A) a/2 f e LP. Here A is the Laplacian and 
(>. \ (-A) al f is defined on the Fourier transform side by \i;\ a f{%). If < a < n this means that 

/ is a function in V which is the Riesz potential of order a of some other function g in LP, 
namely / = c n l/\x\"~ a * g. If a is integer, then W a ' p is the usual space of those functions 
in LP such that all distributional derivatives up to order a are in LP. 

To convey a feeling about the nature of our condition we first discuss the case a = 1. 
Consider the square function 



><: s(/) 2 (x) = 

sJ ■ Jo 



Ib(x, t) - fix) 



2 dt 

— , X 6 



Here / is a locally integrable function on R" and denotes the mean of / on the open 
ball with center x and radius t. One should think of ^ Bu,) f{x) as a quotient of increments 
of / at the point x. Our characterization of W l ' p reads as follows. 



Theorem 1. Ifl < p < oo, then the following are equivalent. 

(1) f€ W hp 

(2) feL p andS(f)eLP. 

If any of the above conditions holds then 

lis(/)ll p -liv/v 



1 



The symbol A B means, as usual, that for some constant C independent of the 
relevant parameters attached to the quantities A and B we have C l B<A<CB. 

Notice that condition (2) in Theorem Q] above is of a metric measure space character, 
because only involves integrals over balls. It can be used to define in any metric measure 
space X a notion of Sobolev space W l,p (X). It is not clear to the authors what are the 
relations of this space with other known notions of Sobolev space in a metric measure 
space, in particular with those of Hajlasz flHll or Shanmugalingam (Si (see also HHKIO . 

The proof of Theorem Q] follows a classical route (see UStrlO . The relevant issue is the 
necessary condition. First, via a Fourier transform estimate we show that 



for good functions /. In a second step, we set up a singular integral operator T with values 
in L 2 (dt/t) such that 



The kernel of T turns out to satisfy Hormander's condition, so that we can appeal to a well 
known result of Benedek, Calderon and Panzone HGRl Theorem 3.4, p. 492] on vector 
valued Calderon-Zygmund Theory to conclude the proof. The major technical difficulty 
occurs in checking Hormander's condition. 

The proof extends without pain to cover orders of smoothness a with < a < 2. The 
square function S (/) has to replaced by 



The result is then that, for < a < 2, / e W a ' p is equivalent to / e If and S a (f) e LP. 
Notice that 



where the barred integral on a set stands for the mean over that set. Stricharzt ( UStrlO used 
long ago the above square function for < a < 1 to characterize W a,p . However the 
emphasis in UStrM was on a larger variant of S a (f) in which the absolute value is inside the 
integral in y in (OQ). In the interval 1 < a < 2 putting the absolute value inside the integral 
destroys the characterization, because then one gives up the symmetry properties of R". 
For instance, S a (j) vanishes if / is a first degree polynomial. 

There are in the literature square functions very close to ® which characterize W a ' p , 
for < a < 2. For example, first differences of / may be replaced by second differences 
and the absolute value may be placed inside the integral ( UStrM and (Stl Chapter V]). The 
drawback with second differences is that they do not make sense in the setting of metric 
measure spaces. See also the paper by Dorronsoro llDll . 

We now proceed to explain the idea for the characterization of W 2 ' p . Take a smooth 
function / and consider its Taylor expansion up to order 2 around x 



P(/)|| 2 = c||V/|| 2 , 



\\T(f)\\i?®«,m m = \\S<j)h- 





(1) 



f(y) = fix) + V/(x) • (y - x) + ^ dPf(x)(y -xf + R, 



(2) 



iei=2 



2 



where R is the remainder and B a multi-index of length 2. Our goal is to devise a square 
function which plays the role of S i(/) (see CQ) for a = 1) with respect to second order 
derivatives. The first remark is that the mean on B(x, t) of the homogeneous polynomial 
of degree 1 in © is zero. Now, the homogeneous Taylor polinomial of degree 2 can be 
written as 

dn f(x) a It 

^(y -xf = H{y -x) + —Af(x) \y - x\ 2 , (3) 
B\ 2n 



for a harmonic homogeneous polynomial H of degree 2. Hence the mean on B(x, t) of the 
homogeneous Taylor polinomial of degree 2 is 



{ 

JB(x, t) 



^-Af(x) \y-x\ 2 dy. 
2n 



This suggests defining 



= f f 

Jo JB(x, t) 



(f(y)-f(x)-±(Af) B(XJ) \y-x\ 2 ) 



dy 



dt 
t 



— , x e 



(4) 



We cannot replace (Af) BixJ ) by Af(x) in the preceding definition, because the mean guar- 
antees a little extra smoothness which one needs in a certain Fourier transform computa- 
tion. Notice that, according to the remarks made before on the mean on the ball B(x, t) 
of the homogeneous Taylor polynomials of degrees 1 and 2, in the expression above for 
S2(f)(x) one may add the missing terms to get the full Taylor polynomial of degree 2, ex- 
cept for the fact that Af(x) should be replaced by (Af) B ( X , t ). Were / smooth enough, one 
could even add the homogeneous Taylor polynomial of degree 3, because it is odd (taking 
x as the origin) and thus its mean on B(x, t) vanishes. This explains why whatever we can 
prove for a = 2 will also extend to the range 2 < a < 4 by considering an appropriate 
square function, which turns out to be 



S «(/)(*)' 



= r f 

Jo JB(xJ) 



(f(y)-f{x)-±Af{x)\y-x\ 2 ) 



dy 



dt 

— , X € 



One should remark that in the range 2 < a < 4 the mean (Af) B ( X ,t) can safely be replaced 
by Af(x) . Here is our second order theorem. 

Theorem 2. If 1 < p < oo, then the following are equivalent. 

(1) f e W 2 -p 

(2) f e LP and there exists a function g e LP such that Siij, g) e LP , where the square 
function S i(f, g) is defined by 



Si(f,g)(x) 2 



r f 

J() JB(x, t) 



(f(y) - fix) - gB(xj) \y - x\ 2 ) 



dy 



dt 

— , x e 



3 



Ifft W 2 ' p then one can take g = A//2n and if (2) holds then necessarily g = A//2n, 



a. e. 



If any of the above conditions holds then 



||S(/,S)|L-||A/|| 



p • 



Notice that condition (2) in Theorem [2] only involves the Euclidean distance on R" and 
integrals with respect to Lebesgue measure. Thus one may define a notion of W 2 ' P (X) on 
any metric measure space X. For more comments on that see section 4. 

Again the special symmetry properties of R" play a key role. For instance, S 2 anni- 
hilates second order polynomials. Theorem [2] has a natural counterpart for smoothness 
indexes a satisfying 2 < a < 4. The result states that a function / e W a ' p if and only if 
/ e LP and there exists a function g & LP such that S a {f,g) e L p , where 



S a (f,g)(x) 2 



r f 

JO JB(x, t) 



{fiy)-f{x)-g{x)\y-x\ 2 ) 



dy 



dt 

— , x e 



Notice that in the range 2 < a < 4 we do not need to replace g(x) by the mean gB( x ,t)- 

Before stating our main result, which covers all orders of smoothness and all p with 
1 < p < 00, we need to discuss a couple of preliminary issues. The first is the ana- 
logue of © for homogeneous polynomials of any even degree. Let P be a homogeneous 
polynomial of degree 2 j. Then P can be written as 

P(x) = H(x) + A j P \x \ 2j , 

where L, = A j (\x \ 2j ) and H satisfies A } H = 0. This follows readily from flSfl 3.1.2, p. 69]. 
Considering the spherical harmonics expansion of P(x) we see that f M=1 H(x)dcr = 0, 

<t being the surface measure on the unit sphere, and thus that [ H(x) dx = 0, t > 0. The 
precise value of Lj, which can be computed easily, will not be needed. 

Our main result involves a square function associated with a positive smoothness index 
a. Let Af be the unique integer such that 2Af < a < IN + 2 . Given locally integrable 
functions /, gi, . . . , g N we set 



S a (f,gl,g2,---,gN)(X) 



r f 

Jo JB(x, t) 



R N (y, x) 



dy 



dt 
t 



X G 



where 



R N (y, x) = f(y) - f(x) - gl (x) \y-x\ 2 g N ^(x) \y - x\ 2{N ~ l) - (g N ) B(x ,t)\y - x\ 2N 

if a = 2N , and 

R N (y, x) = f(y) - f(x) - 8l (x) \y-x\ 2 g N ^(x) \y - x\ 2(N ~ l) - g N (x) \y - x\ 2N 

if2N < a < 2N + 2. 



4 



Theorem 3. Given a > let N be the unique integer such that 2N < a < 2N + 2. If 
1 < p < oo, then the following are equivalent. 

(1) fe W a ' p 

(2) f&LP and there exist functions g j€L p ', 1< j<N, such that S a {f,g\,g 2 , ■ ■ ■ ,gN)€L p . 
If f e W a ' p then one can take gj = A j f/Lj and if (2) holds then necessarily gj = 

A j f/Lj a. e. 

If any of the above conditions holds then 

\\S a (f,gu...,g N )\\ P ~\\(-AT /2 f\\ P . 

Again condition (2) in Theorem [3] only involves the Euclidean distance on R' ! and 
integrals with respect to Lebesgue measure. Thus one may define a notion of W a ' p (X) for 
any positive a and any 1 < p < oo on any metric measure space X. For previous notions 
of higher order Sobolev spaces on metric measure spaces see ULLWL See section 4 for 
more on that. 

The proof of Theorem |3] proceeds along the lines sketched before for a = 1. First we 
use a Fourier transform computation to obtain the relation 

\\S a (f, Af/L u A N f/L N )\\ 2 = c ||(-Ar /2 /ll2 . 

Then we introduce a singular integral operator with values in L 2 (dt/t 2a+i ) and we check 
that its kernel satisfies Hormander's condition. 

The paper is organized as follows. In sections 1, 2 and 3 we prove respectively The- 
orems [TJ [2] and In this way readers interested only in first order Sobolev spaces may 
concentrate in section 1 . Those readers interested in the main idea about jumping to orders 
of smoothness 2 and higher may read section 2. Section 3 is reserved to those interested 
in the full result. In any case the technical details for the proof of Theorem Q] are somehow 
different of those for orders of smoothness 2 and higher. The reason is that Hormander's 
condition involves essentially taking one derivative of the kernel and is precisely the ker- 
nel associated to the first order of smoothness that has minimal differentiability. 

Our notation and terminology are standard. For instance, we shall adopt the usual 
convention of denoting by C a constant independent of the relevant variables under con- 
sideration and not necessarily the same at each occurrence. 

If / has derivatives of order M for some non-negative integer M, then V M f = {d^f)\p\ =M 
is the vector with components the partial derivatives of order M of / and \V M f\ its Eu- 
clidean norm. 

The Zygmund class on R" consists of those continuous functions / such that, for some 
constant C, 

I/O + h) + f(x -h)- 2f(x)\ < C \h\, x,he R" . 

The basic example of a function in the Zygmund class which is not Lipschitz is fix) = 
\x\ log |x|, x eW\ 

The Schwartz class consists of those infinitely differentiable functions on R" whose 
partial derivatives of any order decrease faster than any polynomial at oo. 
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After the first draft of the paper was made public Professor Wheeden brought to our 
attention his articles HWll and HW2L In HWll a general result is proven which, in partic- 
ular, contains Theorem 3 for < a < 2. In HW2II Sobolev spaces with respect to special 
homogeneous spaces are considered for < a < 1 . 



2 Proof of Theorem U 

The difficult part is the necessity of condition (2) and we start with this. 
As a first step we show that 

H5 , i(/)|| 2 = c||V/|| 2 (5) 

for a dimensional constant c. Set 

x{x) = \mT)\ XBm){x) 

and 

1 x 

Xt(x) = — . 

so that, by Plancherel, 

S l (f)(x) 2 dx = \(f*Xt)(x)-f(x)\ 2 dx- 
Jo Jr" ' 



r r mz)-i\ 2 \m z d^. 

Jo Jr» * 



Since % is radial, x(€) = F(\g\) for a certain function F defined on [0, oo). Exchange 
the integration in d£ and dt in the last integral above and make the change of variables 
r = t\%\. Then 

f S x (f)(x) 2 dx = c [ f |F(r) - 1 1 2 ^ |/Wl£| 2 
Jr" Jr» Jo t 



= c | |F(r)-l| 2 ^ IIV/H 2 



and © is reduced to showing that 



f 

Jo 



|F(r)-l| 2 -— <oo. (6) 

Set 5 = 5(0, 1) and g, = (1, 0, . . . ,0) e R". Then 

F(t) =x(te 1 ) = \ expi-ix^dx 
Jb 

1 - ix\t - 2 X i ?2 + " " " 



Jr 



= 1 - ^ -j' x 2 dxt 2 + 
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which yields 

F(t) - 1 = 0(t 2 ), as t -> 

and shows the convergence of © at 0. 

Since F(|£|) = ,£(£) is the Fourier transform of an integrable function, F(t) is a 
bounded function and so the integral © is clearly convergent at oo. 

We are left with the case of a general p between 1 and oo. If / e W l ' p , then / = 
g*l/|jt|" _1 for some g e LP (with l/|jt| n_1 replaced by log |x| for n = 1). Set/(jc) = 1 /|jc|" 1 . 
Then 

~ /(*) = (f*Xt)(x) - fix) = (g * K t ){x) , 

where 

K t (x) = (I * X t)(x) - I(x) = f I(y) dy - I(x) . (7) 

JB(x, t) 

If we let T(g)(x) = (g * K,){x), x e W, then one can rewrite S \(f)(x) as 

i 

J o \(g*K t )(x)\ 2 -J =\\Tg(x)\\ L 2 (dm . 

Then © translates into 

f \\Tg(x)\\* ) dx = c\\g\\l, 

Jr" 

and we conclude that T is an operator mapping isometrically L 2 (R.") into L 2 (K", L 2 (dt/ 1 3 )). 
If the kernel K t (x) of T satisfies Hormander's condition 

f \\K t (x -y) - K t (x)\\ L 2 (dt/t3) dx<C, ye R" 
J\4>2b>\ 

then a well known result of Benedek, Calderon and Panzone on vector valued singular 
integrals (see HGRl Theorem 3.4, p. 492]) yields the L p estimate 

§jT g mi Hdm dx<c p \\g\\i, 

which can be rewritten as 

I|5i(/)IIp<C p ||V/|| p . 

The reverse inequality follows from polarization from © by a well known argument 
( HGR[ p. 507]) and so the proof of the necessary condition is complete. We are going to 
prove the following stronger version of Hormander's condition 

\\K t (x-y)-K t (x)\y idm <C-^, yeR", (8) 

\x\ 

for almost all x satisfying |x| > 2\y\. 
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To prove © we deal separately with three intervals in the variable t. 

Interval 1: t < y. From the definition of K, in © we obtain 

Vtf,(x) = (W ***)(*) -V/(x). 

Notice that, in the distributions sense, the gradient of / is a constant times the vector 
valued Riesz transform, namely 

x 



VI = -(n - l)p.v. 



x \n+\ ■ 



If |jc| > 2\y\, then the segment [x - y, x] does not intersect the ball 5(0, \x |/2) and thus 

\K t (x-y)-K t (x)\<\y\ sup \VK t (z)\. (9) 

zelx-yy] 

If t < \x\/3 and z e [x - y,y], then B(z, t) c R" \ 5(0, \x |/6), and hence 

VK t (z) = f (V/(w) - V/(z)) . (10) 
Taylor's formula up to order 2 for V/(w) around z yields 



w-z 



2 



V/(w) = V/(z) + V7(z)(w - z) + j-) , 

where V 2 /(z)(w - z) is the result of applying the matrix V 2 /(z) to the vector w - z. The 
mean value of V 2 /(z)(w - z) on 5(z, t) is zero, by antisymmetry, and thus, by (fTOl) . 

|VTOI<C 1 



|jc|" +2 
and so, by © 



r 



\K t (x-y)-K t (x)\<C\y\- 

\.x\ 



Integrating in t we finally get 

f w/3 7 dtV m ( r w/3 



tdt\ =C 



\y\ 



x\ n+l 



Interval 2: |x|/3 < f < 2|x|. The function / is continuously differentiable on 
R' ! \ S t , S t = {x : \x\ = t}, because its distributional gradient is given by / * Vxt an d each 
component of Vxt is a Radon measure supported on S ,. The gradient of / * Xt is given at 
each point x e R" \ S , by the principal value integral 



p.v,(VI *Xt)(x) = ~(n - \)p.v. \ -—j 

jB(x,t) \y i 



which exists for all such x. The difficulty in the interval under consideration is that it may 
happen that |x| = t and then the gradient of / * Xt has a singularity at such an x. We need 
the following estimate. 
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Lemma 1. 



p.v. I 

JB{x,t) 



\y\ n+1 



dy 



< C log 



\x\ + t 
\x\-t\'- 



x e 



Proof. Assume without loss of generality that x = (x\, 0, . . . , 0). The coordinates yj, j ± 
1, change sign under reflection around the y\ axes. Hence 



Now, if \x | < t, 



If \x | > t, 



>.v. I -^-r dy = 0, 1 < /' < n . 



JB(x,t) 



y^ 



\y\ 



n+ 1 



<c - 



.r-W) M 



(1+1 



dy 



Clog 



t - \x\ 



p.v. I 

JB(x, t) 



Vl 



Ivl 



«+ 1 



dy 



L( X , t) \y\ n+l d 

: r^'d 1 = 

J\x\-t t 



Clog 



|x| + t 
x\ — t 



Assume without loss of generality that y = (y l5 0, . . . , 0). The distributional gradient 
of/ is 

-in - \)py.-j-r *Xt, 



\y\ 

which is in L 2 . Then / *Xt £ W 1,2 and consequently is absolutely continuous on almost all 
lines parallel to the first axes. Therefore 



K t (x-y)-K t (x) 



Jo 



VK,(x -ry)-ydr 



for almost all x and 



dr . 



Hence 

>2|x| 



/ r2\x\ 
\J\x\/3 



dtV 



\K t {x-y)-K t {x)\ l -\ <C 



■|n+l 



Jw/3 \Jo 



1 +log 



|jc - ry| + ? 
I* - ry\ - t\ 



where the last identity is a definition of D. Applying Schwarz to the inner integral in D 
and then changing the order of integration we get 



D 2 < 



Jo WW/3 

For each r make the change of variables 



1 +log 



|x - ry| + t \ dt 
\\x-Ty\-t\) t 



dr . 



s = 



\x - ry| 



to conclude that 



D 2 < 



2/9 



1 +log 



1 + s \ ds 



l-s\ 



Interval 3: 2\x \ < t. For each z in the segment [x - y,y] we have 5(0, t/4) c B(z, t). 
Then, by (0, 



l^feOI Jfifz.rt |W|" + 



W , Z 

■ dw - 



f)\B(0,f/4) I 14 ' 



n+1 



kl 



and so 



Hence, owing to ©, 



and thus 



\VK t (z)\<C—, ze[x-y,y\. 
\x\ 



\K t (x-y)-K t (x)\<C^l 

\x\" 



f °° 9 ^A 5 Lvl / C°° dt\' 2 \y\ 

\K t (x-y)-K t (x)\ 2 -\ <C^f T = C ^T> 

which completes the proof of the strengthened form of Hormander's condition ©. 

We turn now to prove that condition (2) in Theorem Q] is sufficient for / e W l,p . Let 
/ e LP satisfy S\(f) e LP. Take an infinitely differentiable function cf> > with compact 
support in 5(0, 1), f (f> = 1 and set <p e (x) = pr^(f), £ > 0. Consider the regularized 
functions f £ = f* 4> e . Then f e is infinitely differentiable and ||V/ e || p < ||/|| p ||V0 e ||i, so that 
f £ e W l ' p . Thus, as we have shown before, 

HV/ 6 ||*||Si(/ 6 )|lp. 
We want now to estimate \\S i(fe)\\ P independently of e. Since 

(/ e W,r) - fe(x) = ((/ *Xt~f)* <Pe) (X) , 
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Minkowsky's integral inequality gives 



S l(fe)(x) = \\(f £ ) B ( X ,t) ~ fe{x)\\ mm < (S !(/) * <f> e )( X ) , 

and so ||V/ e || < C ||S i(/)|| p , e > 0. For an appropriate sequence ej — > the sequences d k f € , 
tend in the weak ★ topology of U to some function gk e LP, 1 < k < n. On the other 
hand, f e — » / in L p as 6 — > and thus c^/e — » dfc/> 1 < < n in the weak topology of 
distributions. Therefore d k f = g k for all k and so / e W l,p . 



3 Proof of Theorem [2] 

The difficult direction is (1) implies (2) and this is the first we tackle. We start by showing 
that if / € W 2 ' 2 then 

l|S 2 (/)|| 2 = c||A/|| 2 (11) 

where the square function S 2 (f) is defined in ©. To apply Plancherel in the x variable it 
is convenient to write the innermost integrand in © as 



f ( 
•f ( 



f( x + h)-f(x) 



(f 

WB(x, t) 



In 



dz\\y- xf dy 



(f 

\ J 5(0,0 



Afix + k) 
2n 



dk\\h\ 2 \dh. 



Applying Plancherel we get, for some dimensional constant c, 



c\\s 2 (fm 



r r f 

Jo Jr"Jb(o,o 



exp (igh) - 1 + 



(f " 

\ J 5(0,0 



exp (igk) dk 



2n 



dh\m\ 2 d^ 



dt 
? 



Make appropriate dilations in the integrals with respect to the variables h and k to bring 
the integrals on 5(0, 1). Then use that the Fourier transform of |g i 1 ^ a(o,i) is a radial 
function, and thus of the form F(|£|) for a certain function F defined on [0, oo). The result 
is 



c\\S 2 (f)\\: 



-If 

Jr« Jo 



Fim-i+^F^^-f \h\ 2 dh 

^ n Jb(o.i) 



CLX o 

^l/(£)l 2 ^. 



The change of variables r = / 1£| yields 



where / is the integral 

'"f 

Jo 



c||S 2 (/)ll2 = /||A/|g 



F(r) - 1 + t 2 F(t)^- £ \h\ 2 dh 
2n Js(o,i) 



dr 

~5 • 



(12) 
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The only task left is to prove that the above integral is finite. Now, as r — » 0, 

F(t) = t exp(ihiT)dh 

Jfi(0,l) 



f ( 

^8(0,1) \ 



1 + ih X T - ]^h 2 \T 2 + • • • 



l-iff /z 2 J/z)t 2 + 0(t 4 ). 

^ Vb(o,i) 



Hence 



F(r) - 1 + t 2 F(t)^- f \h\ 2 dh 
2^ Js(o,i) 



LiS h^dh + LS \h\ 2 dh)r 2 + 0(T 4 ) = 0(T 4 ), 

\ Jfl(0,l) ^ n Jb(0,1) ) 



because clearly f \h\ 2 dh = n f h 2 dh. Therefore the integral (PT21) is convergent 
at t = 0. 

To deal with the case t -> oowe recall that F can be expressed in terms of Bessel 
functions. Concretely, one has ( llGrl Appendix B.5, p. 429]) 

P,iP(r) = ^. 
The asymptotic behavior of /„/2(t) gives the inequality 

\F(r)\<C^<ci 

T — T 

which shows that the integral (PT21) is convergent at oo. 

We turn our attention to the case 1 < p < oo. Let I 2 (x) stand for the kernel defined on 
the Fourier transform side by 

1 

l^ 2 

In other words, I 2 is minus the standard fundamental solution of the Laplacian. Thus 
h(x) = c n l/\x\ n ~ 2 if n > 3, I 2 (x) = log |x| if n = 2 and I 2 {x) = -\ \x\ if n = 1. Given 
any / e W 2 ' p there exists g e LP such that f = I 2 * g (indeed, g = -Af). We claim that 
there exists a singular integral operator T(g) taking values in L 2 (dt/t 5 ) such that 

S 2 (f)(x) = \\T(g)(x)\y (dt/fi) . (13) 

Set 

1 

x{x) ~ \mv)\ XBm){x) 



= 77£ ■ 
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and 

1 x 
Xt(x) = -X(-)- 



Then, letting M = f \z\ 2 dz, 



f ( 



B(0, 1) 

f(y) - fix) - ^(A/W) \y - x\ 2 \ dy = ((h * Xt -h-^ t 2 X t) * 8)(x) 



= (K t *g)(x), 

where 

K t (x) = (I 2 * Xt )(x) - I 2 (x) - ^ t 2 xt(x) . 

Setting T(g)(x) = (K t * g)(x) we get from the definition of 5 2 (/) in ©. Then (HB 
translates into 

\\T g(x)\\ 2 mm dx = c\\g\\l, 

and we conclude that T is an operator mapping isometrically L 2 (W) into L 2 (W, L 2 (dt/t 5 )), 
modulo the constant c . If the kernel K t (x) of T satisfies Hormander's condition 

f HZ"^ - y) - K t (x)\\ L 2 (dtm dx<C, ye R", 
^W>2|yl 

then a well known result of Benedek, Calderon and Panzone on vector valued singular 



integrals (see [GR, Theorem 3.4, p. 492]) yields the LP estimate 



\\Tgm>dx<C P \\g%- 



which can be rewritten as 

IIW)||p<C p ||A./V 

The reverse inequality follows from polarization from (fTTT) by a well known duality argu- 
ment (" HGR1 p. 507]) and so the proof of the necessary condition is complete. 

We are going to prove the following stronger version of Hormander's condition 

\\K t (x - y) - K t {x)\\ L 2 {dm < C \x\ > 2\y\ . 

For this we deal separately with the kernels H,(x) = (I 2 * Xt){x) - h{x) and t 2 Xt{x)- For 
t 2 Xt(x) we first remark that the quantity \Xt(x-y) ~Xt(x)\ is non-zero only if \x-y\ < t < \x\ 
or \x\ < t < \x—y\, in which cases takes the value 1 /c n t", c n = \B(0, 1)|. On the other hand, 
if \x\ > 2\y\ then each z in the segment joining x and x-y satisfies |z| > \x\/2. Assume that 
|x - y\ < \x\ (the case \x\ < \x - y\ is similar). Then 

i 1 



3 , , dt\' 2 i r w dt V 



., -, <_c M" 2 



\x-y\ 2n \x\ 2n j |x|" +1/2 
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We check now that H, satisfies the stronger form of Hormander's condition. If t < 
\x |/2, then the origin does not belong to the ball B(x - y, t) nor to the ball B(x, t). Since 
I 2 is harmonic off the origin, the mean of I 2 on these balls is the value of I 2 at the center. 
Therefore H t {x - y) - H t (x) = in this case. 

If f > \x\/2, then 

\H t (x-y)-H t (x)\<\y\ sup \VH t (z)\ < C . 
The last inequality follows from 

V// ( (z)=f VI 2 (w)dw-VI 2 (z), 
|V/ 2 (z)| <Cl/\z\"~ l ^Cl/lxl^and 

if V/ 2 (w)Jw|< f JL^c* 

JB(z,t) JB(z,t) \ w \ Kl 

We turn now to prove that condition (2) in Theorem [2] is sufficient for / e W 2 ' p . 
Let / and g in V satisfy S 2 (f,g) e ZA Take an infinitely differentiable function <p > 
with compact support in 5(0, 1), § <p = 1 and set <p e {x) = pr^(f), e > 0. Consider the 
regularized functions f e = f * <Pe an d g £ = g * 4> £ - Then / e is infinitely differentiable and 
l|A/ e || p < ll/llpllA^Hi, so that f e G W 2 A Recalling that M = £ (01) kl 2 *, we get, by 
Minkowsky's integral inequality, 



Therefore 



S 2 (fe,ge)(x) = \\(f £ *XtXx)-Mx)-(g e *XtXx)M*t*\\ lHdt/fi) 

= II ((/ *xt)-f-ig *xd m 2 f- * <f> e ) (x)\y (dt/fi) 

<(S 2 (f,g)*<f> e )(x). 
Now we want to compare (1 /2n)A/ e and g e . Define 

(r°° 1 2 dt^^ 2 

I M 2 — (A/ e **,)(*) 



Then 

DAx)<S 2 (f £ )(x) + S 2 (f„g £ )(x) 

<S 2 (fMx) + (S 2 (f,g)(x)* ( f> £ )(x), 
and thus Z) £ (jc) is an LP function. In particular D £ (x) < 00, for almost all x e W. Hence 
|(l/2n)A/ e (x) -g £ (x)\ = lim|(l/2n)(A/ e * Xt )(x) - (g e * Xt )(x)\ = 0, 

r->0 

for almost all 1 G R", and so (l/2«)A/ e — > g in Z/ as e — > 0. Since f e — > / in Z/ as 6 — » 0, 
then A/ e — » A/ in the weak topology of distributions. Therefore (l/2n)A/ = g and the 
proof is complete. 
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4 Proof of Theorem ffl 



The difficult direction in Theorem[3]is to show that condition (2) is necessary for / 6 W a ' p . 
The proof follows the pattern already described in the preceding sections. One introduces 
an operator T taking values in L 2 (dt/t 2a+l ) and shows via a Fourier transform estimate that 
T sends L 2 (R") into L 2 (R.",L 2 (dt/t 2a+l )) isometrically (modulo a multiplicative constant). 
The second step consists in showing that its kernel satisfies Hormander's condition, after 
which one appeals to a well known result of Benedek, Calderon and Panzone on vector 
valued singular integrals to finish the proof. 



4.1 The fundamental solution of (-A) a/2 

Let I a be the fundamental solution of (—A) a/2 , that is, I a is a function such that I a (£) = \£\~ a 
and is normalized prescribing some behavior at oo. It is crucial for our proof to have an 
explicit expression for I a . The result is as follows (see HACL1 or HMOPVl p. 3699]). 
If a is not integer then 

Ia(x) = C a , n \x\°-'\ (14) 

for some constant c a>n depending only on a and n. 

The same formula works if a is an even integer and the dimension is odd or if a is an 
odd integer and the dimension is even. 

The remaining cases, that is, a and n are even integers or a and n are odd integers are 
special in some cases. If a < n formula (PT4l) still holds, but if a is of the form n + 2N, for 
some non-negative integer N, then 

4 (*) = c a , n \x\ a -" (A + B log \x |), x 6 R" , 

where c ain , A and B are constants depending on a and n, and 5^0. Thus in this cases 
(and only in this cases) a logarithmic factor is present. For instance, if a = n, then 
I a (x) = B log \x |. If n = 1 and a = 2, then I 2 (x) = -(1/2) |x | and there is no logarithmic 
factor. 



4.2 The case p = 2 

Given a positive real number a let iV be the unique integer satisfying 2N < a < IN + 2. 
Define the square function associated with a by 



= r f 

Jo JB(x, t) 



P N(y, x) 
t a 



dy 



dt 

— , X G 



(15) 



where p N (y, x) is 



f(y)-f(x)-^Af(x) \y-x\ 
2n 



A"-7(*)|y-jc|- 



^-(A"f) B ( X Jy-x[ 
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if a = 2N , and 

f(y) ~ fix) - -U/(x) \y - x\ 2 J-A^VCx) \y - x\ 2 ^ r) - -}-A N f(x) \y - x\ 2N 

if 2N < a < 2N + 2 . 

Recall that Lj = A j (\x \ 2j ) and that the role which the L ; play in Taylor expansions was 
discussed just before the statement of Theorem [3] in the introduction. 
In this subsection we prove that 

l|5 ff (/)|| 2 = c||(-Ar /2 (/)ll2. (16) 

We first consider the case a = IN and then we indicate how to proceed in the (simpler) 
case IN < a < IN + 2 . Our plan is to integrate in x in (fT5T) . interchange the integration 
in x and t and then apply Plancherel in x. Before we remark that by making the change of 
variables y = x + th we transform integrals on B(x, t) in integrals on 5(0, 1) and we get 

£ p N (y,x)dy=£ fix + th) dh-Y t 2} £ \h\ 2j dh 

- £ A N f(x + th) dh t 2N £ \h\ 2N dh . 

JB(0, 1) Jb(0, 1) 

Now apply Plancherel in x, as explained before, and make the change of variables r = t 
where £ is the variable in the frequency side. We obtain 

l|S«(y)lll = c/||(-A) o/2 /lll, (IV) 

where 



1= F(T)-V(-l)¥^-(-l)V^F(T) 
Jo ^ L N 

Here F is a function defined on [0, oo) such that F(\g\) gives the Fourier transform of 
the radial function * ^ 8(0,1) at tne point and we have introduced the notation Mj = 
f j \h\ 2j dh. We have to show that the integral / is finite. 

Using the series expansion of the exponential we see that, as r — > 0, 

F(t) = T exp(ihir)dh 

Jfi(0,l) 

= 1 + • • • + (-1)^ T 2N -±- £ hfdh + 0(T 2N+2 ) . 

(zyvj! Jb(o,i) 

We need to compare f B(Q hf 'dh with f m 1} \h\ 2N dh. The linear functionals P -» A 2j (P) 
and P — > ^£ (Q P, defined on the space #27 of homogeneous polynomials of degree 2 j, 



2 

dr 



T 2a+1 
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have the same kernel. This follows from the discussion before the statement of Theorem^] 
in the introduction. Therefore, for some constant c, 



A 2j (P) = c £ 



P, Pt H 2j . 

B(0,1) 



Taking P{x) = \x \ 2j we get L, = c ■£.. „ \x\ 2J = dx, and taking P{x) = x 2j we get 
=c fw,»*? dx - Hence 



-7')! Jb(o,i) Jb(0,1) Lj 



(2y)! 

and thus, owing to the definition of / and the fact that F(r) = 1 + 0(t 2 ), as r — » 0, 

Then 7 is convergent at because a < IN + 2 (indeed, now a = IN , but this part of the 
argument works for the full range IN < a < IN + 2 .) 

We turn to the case r — > oo. Notice that the only difficulty is the last term in the 
integrand of I, because 



X 



T 



: dr 



T-2a+ 
1 T 



r <CX3, < j < N ~ 1 , 



provided 2/V < a. To deal with the term 



t 2N F(t)\ 2 4-T (18) 



we only need to recall that F can be expressed in terms of Bessel functions. Concretely, 
one has (EjJ Appendix B.5, p. 429]) 

|S(0,l)|F(r) = ^. 

The asymptotic behaviour of /„/2(t) gives the inequality, as r — > oo, 

|F(r)|<c4r^ci 

which shows that the integral (PT8l is finite provided 2iV - 1 < a, which is the case because 
a = IN . Observe that the argument works for IN < a < IN + 2 provided pNiy, x) is 
defined replacing A N f(x) by (A N f) B ( X ,t)- 

In the case 2N < a < IN + 2 we argue similarly. After applying Plancherel we obtain 
(fTTT) where now I is 



i= f(t) - Yi-iy^j-L 



2 

dr 



T 2»+l 
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The proof that / is finite as r — » is exactly as before. As r — > oo the situation now is 
simpler because the worst term in / is the last one, namely, 



Ji 



T 4/V ^ T 



T 2a+1 



which is finite because 2N < a. 



4.3 A vector valued operator and its kernel 

Given / e W a ' p , there exists a function g e If such that / = I a *g. Indeed, g = (-A) a/2 (f). 
Then 

\ p N (y, x) dy = (K t * g)(x) , 

JB(x, t) 

where the kernel K t (x) is 
I 



= f I a (y) - V - A'/ a (*) |y - x| 2j - — (A\) fi( , d |y - x| 



2N 



dy (19) 



if a = 2jV , and 



K,(x) 



r ( N i 1 

f ^(y)-Y! y A j I a (x)\y-x\ 2j 
JB(x,t){ j=0 h ) 



dy 



(20) 



if IN < a < IN + 2 . 

Hence the square function associated with the smoothness index a is 

S a (f) 2 (x) = \(K t * g)(x)\ 2 — , 
Jo ' 

Define an operator T acting on functions / e L 2 (R") by 

Tg(x) = (K t * g)(x), xel". 

The identity (HSJ in subsection 4.2 says that T takes values in L 2 (R n ,L 2 (dt/t 2a+l )) and, 
more precisely, that 

\\Tg(x)\\ 2 L2{dt/fla+l) dx = \\S a (f)\\l = c \\g\\l . 

Therefore T is an operator mapping isometrically (modulo a multiplicative constant) 
L 2 (R") into L 2 (R", L 2 (dt/t 2a+i )) and we have an explicit expression for its kernel. If we 
can prove that K t (x) satisfies Hormander's condition 

f \\K t (x - y) - K t (x)\\ L 2 (dt/t 2„ +1) dx <C, yeW 1 , 

J\x\>2\y\ 
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then the proof is finished by appealing to a well known result of Benedek, Calderon and 
Panzone (EH Theorem 3.4, p. 492]; see also (GRl p. 507]). In fact, we will show the 
following stronger version of Hormander's condition 

\\K t (x -y)- K t (x)\\ L 2 (dt/t 2 a+1) < C \x\ > 2\y\ , (21) 

\x\ 

for some y > depending on a and n. 

The proof of (I2TT) is lengthy. In the next subsection we will consider the case of small 
"increments" in t, namely t < |jc|/3. 

4.4 Hormander's condition: t < \x\/2> 

We distinguish two cases: 2N < a < 2N + 2 and a = IN. Assume first that IN < a < 
2N + 2. Set 

x{x) = whs\ xmi)(x) 

and 

1 x 
Xt(x) = -Xi~) ■ 

To compute the gradient of K, we remark that 

K t (x) = (4 *Xt)(x) -YjT^t AJI ^ x) ' 
j=o L i 

where M ; - = f m \z\ 2j dz. Thus 



r ( N i "1 

Vtf f (*) = f V4(v) - V - A j (VI a )(x) \y - x\ 2j 

JB(x,t){ j=0 L j ) 



dy. 



Let P m (F, x) stand for the Taylor polynomial of degree m of the function F around the 
point x. Therefore 

VK t (x) = f (VI a (y) - / WV4, x)(y)) dy , 

JB(x, t) 

because the terms which have been added have zero integral on the ball B(x, t), either 
because they are Taylor homogeneous polynomials of VI a of odd degree or because they 
are the "zero integral part" of a Taylor homogeneous polynomial of VI a of even degree 
(see the discussion before the statement of Theorem 3 in the introduction). Given x and y 
such that \x\ > 2\y\, apply the formula above to estimate VK t (z) for z in the segment 
from x - y toy. The standard estimate for the Taylor remainder gives 

\VK t (z)\ < t 2N+2 sup |V 2W+3 Uw)|. 

weB(z, t) 
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Notice that if z £ [x - y,y], w e B(z, t) and t < |x|/3, then \w\ > \x\/6. Now, one has to 
observe that 

iv^^i^ciwr 1 - 2 ^ 3 , 

owing to the fact that logarithmic factors do not appear because the exponent a - n - IN - 
3 < — n - 1 is negative. By the mean value Theorem we then get 



\K t (x-y)-K t (x)\<\y\ sup \VK t (z)\ < C \y\t 2N+2 \x\ a -"- 2N ^ . 

Z€[x-y,y\ 



Since 



we obtain 



I 



w/3 dt \ 1/2 

f 2(2N+2) Ul | _ q \x\ 2N+2 ~ a 
*2or+l 
* 



' r yi 7 dt \ 1 \y\ 



which is the stronger form of Hormander's condition (12T1) with 7= 1 in the domain t< \x\/3. 

Let us consider now the case a = IN. Since A n I 2 n is the Dirac delta at 0, A n I 2 n(x) = 0. 
Hence K,(x) = K ( t l \x) - K (2 \x), where is given by (|20]) with a replaced by IN and 

K?\x) = f -L( (A^) B(;c , f) b - x\ 2N dy = ^ t 2N (A N I 2N ) B(XJ) . 

JB(x,t) L n 

The kernel K\ } is estimated exactly as in the first case by just setting a = IN. The kernel 
requires a different argument. 

Since A n I 2 n is the Dirac delta at the origin, K {2) is a constant multiple of t 2N Xt- We 
show now that this kernel satisfies the strong form of Hormander's condition. The quantity 
\Xt{x-y) ~Xt(x)\ is non-zero only if \x-y\ < t < \x\ or \x\ < t < \x-y\, in which cases takes 
the value 1 /c n f, c n = \B(0, 1)|. On the other hand, if |jc| > 2\y\ then each z in the segment 
joining x and x - y satisfies |z| > \x\/2. Assume that \x - y\ < \x\ (the case \x\ < \x - y\ is 
similar). Then 



1 



LV <c^l 

\x-y\ 2n \x\ 2n ) ~ Ix^ 1 ' 2 ' 



which is (ED with y = 1/2. 



4.5 Hormander's condition: t > \x\/3 

We distinguish three cases: a<n + l,a = n + l and a > n + 1. 
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If a < n + 1, then all terms in the expressions (TT91 ) and (12T)1) defining K t satisfy 
Hormander's condition in the domain t > |x|/3. If a = IN the last term in (PT9l) is of 
the form 



-f 



1 



M 



— \y - xY N dy = --^ t 2N (A N I a ) B(xJ) = Ct^xlx) , 

L N 



B(x,t) Ln 



which has been dealt with in the previous subsection. Let us consider the terms of the 
form t 2j A j I a (x), j > 0. One has the gradient estimate 



\t 2j VA j I a (x)\ < C t 2i \x\ a - n - 2j - 1 , 



(22) 



because no logarithmic factors appear, the reason being that the exponent a — n — 2j— \ < 
a - (n+ 1) is negative. Since 



WW/3 t 2a+l ) 



C \x 



2j-a 



we get Hormander's condition with y = 1 in the domain t > \x |/3. 

It remains to take care of the first term £ I a (y) dy in (|T9l ). We have the following 
obvious estimate for its gradient 



f V4(y)J<cf \yr n - l dy. 

JB(x,t) JB(x,t) 



Notice that there are no logarithmic factors precisely because a < n + 1 . The integrand 
in the last integral is locally integrable if and only if a > 1. Assume for the moment that 
1 < a < n + 1 . Then 

4—n+a—l 



j 

JB(x, t) 



VI a (y)dy 



< Ct 



Since 



if 

WM/3 



,2(a-n-l) 



dt 

t2a+\ 



1/2 



C \x\ 



-n-l 



we get Hormander's condition with y = 1 in the domain t > \x |/3. The case a = 1 has 
been treated in section 1, so we can assume that < a < 1. We need the following well 
known and easily proved inequality 

Lemma. Let E be a measurable subset ofW and < J3 < n. Then 

dz 



Ie \z 



n-p 



< C \Ef" , 



where \E\ is the Lebesgue measure ofE. 
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Denoting by D the symmetric difference between B(x, t) and B(x-y, t), we obtain, by 
the Lemma, 



f I a (y)dy-f I a (y)dy\<Cr" f JjjjL 

JB(x-y,t) JB{x,t) JD U\ 



Since 



if 

WM/3 



.2(a—n—alri) 



<Ct n (t"- 1 \y\) a/ " = Cf- n ~ a,n \y\ aln . 
dt 



t 2a+\ 



1/2 



C \x\ 



-n-a/n 



we get Hormander's condition with y = a/n in the domain t >\x |/3. 

We tackle now the case a = n + 1. Since a and n are integers with different parity no 
logarithmic factor will appear in I a . Thus I a (x) = C\x\. The proof above shows that the 
terms t 2j AU a (x) appearing in the expression ( fT9l ) of the kernel K, still satisfy Hormander's 
condition for j > 1 . The remaining term is 



f ' 

JB(x, t) 



(I a (y)-I a (x))dy 



and its gradient is estimated by remarking that the function \x | satisfies a Lipschitz condi- 
tion. We obtain 



But clearly 



f 1 

JB(x, t) 



ff *) 

\J\x\/3 t 2a+1 J 



(V400 - V4(x)) 



< c 



1/2 



= c ur = c 1*1 



-n-l 



which completes the argument. 

We turn our attention to the case a > n + l.lf a = 2N the part of K t which has to be 
estimated is 



r N ~ l i 

H,(x) = f I a (y) - V — A J 4W |y - x 



dy . 



Since A 2A 7 ff is the Dirac delta at the origin and x ^ we have A 2N I a (x) = and so H,(x) 
coincides with the expression (1201) for K t (x) in the case 2// < a < 2N + 2 . We are then 
going to deal with this kernel in the full range 2N < a < IN + 2. Let M be the unique 
positive integer M such that -1 < a - n - 2M < 1. We split K, into two terms according 
to M, that is, K, = - K?\ where 



*?>(*) 



r ( M-l \ 

T 4W-y r A^(x)|y-x| 2 ^ 



dy 



and 



* f (2) « 



f VlA^ a (x)|y-x| 
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The estimate of each of the terms in K ( t 2) is performed as we did for the case a < n + 1 
(if a = 2N, then j < N - 1). The gradient estimate is exactly (1221) . Now no logarithmic 
factors appear because the exponent satisfies a - n - 2j — \ < a - n - 2M - 1 < 0. The 
rest is as before. 

To estimate we distinguish three cases: -1 < a-n-2M < 0, < a-n-2M < 1 
and a - n - 2M = 0. In the first case we write the gradient of K$' as 



M-i -, ^ 

dy 

B(x,t){ *p^H 



VK ( t l \x) = ( W a (y) - V ]- A j (VI a )(x) \y - x) 2j 
Jb(xj) L i 



= f (V/ a (y) - />2M-2(V4, x)(y) ) , 

JB(x, t) 

where Pjm-i is the Taylor polynomial of degree 2M - 2 of V/ a around the point x. As 
before, the added terms have zero integral on B(x, t) either because they are homogeneous 
Taylor polynomials of odd degree or the "zero integral part" of homogeneous Taylor poly- 
nomials of even degree. Now fix y in B(x, t) but not in the half line issuing from x and 
passing through the origin. Define a function g on the interval [0, 1] by 

g(T) = V40 + T (y - x)) - P 2 m- 2 (V4, x)(x + r(y- x)), < r < 1 . 

Then g e C°°[0, 1] because the segment with endpoints x and y omits the origin. Since 
^(0) = 0, < j < 2M - 2, 

2 ^ 2 ^(0) 
V/«(y) - F 2 m- 2 (V4, x)(y) = g(l) - £ 



;'=o 7 



Jo 



(2M-2)! * W ' 



by the integral form of Taylor's remainder. The obvious estimate for the derivative of g of 
order 2M - 1 is 



| g 2M-!) (T)| < |V 2M-! ^ + T(y _ my _ X| 2M-1 < c 



x + r(y - ^)|»-(«-2M) 



Since we are in the first case, a is not integer and thus no logarithmic factor exists. More- 
over < n — (a - 2M) < 1, which implies that and that \/\z \ n ~ ia ~ 2M ^ is locally integrable 
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in any dimension. Therefore 

.1 



|V^ f (1) (jc)| < C f lt 2M - 1 -" f dy s , ou ] dr, 

' Jo I JB( xA \x + T(y-x)r^-W) 

c ZM-l-n f l l f \ dT 

Jo \J B (,,r)kr (ff - 2M) /r" 
< c? 2M-i-„ f\ trr -2Mdr 

Jo T" 

= ? «-«-i f 1 tL = C f 

Jo r»-(- 2M ) 



t a—n—\ 



Since 



(f 

Vw, 



dt \ 1/2 



2«+l 



W/3 

we get Hormander's condition with y = 1 in the domain ? > \x |/3. 

Let us consider the second case: < a-n-2M < 1. This time we express the gradient 
of£ f (1) by means of a Taylor polynomial of degree 2M - 1 : 



r I i 

V^ (1) (x) = f V/ a (y) - V - A'(V/ a )(*) \y - x\ 2j 

JB(x,f)\ j^O j 



dy 



= f (V4(y) - P 2M -i(V4, x)(y) ) . 

Using again the integral form of the Taylor remainder of the function g, with P2M-2 re- 
placed by P 2 m-i, we obtain 

-I'd 

JO \JB(x,tT) 



= Ct 2M ~ 



dz \ dr 

^ | z |n-(a-2M-l) j r « 



< 



Cr 2M ~" r\ tT ) a - 2U - 1 — 

Jo T " 



1 dr_ 

T «-(a-2M-l) 



— ___ = c^"- 1 



from which we get the desired estimate as before. 

We turn now to the last case left, a = n + 2M, with M a positive integer. In this case 



I a (x) = C\x\ 2M (A + B\og\x\), xeR", 5*0, 
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where A, B and C are constants depending on n and M. We also have 
tf~%(x) = C |x| 2 (A t + 5! log |jc |), x 6 R" 

and 

VA M ~ l I a (x) = Cx(A 2 + B 2 log |jc |), x e R n . 
In particular VA^" 1 /^ is in the Zygmund class on R". We have 

r ( M-2 ^ 

VK ( t l \x) = \ V/ B (y) - V - A\VI a )(x) \y-x\ 2j - A M ~ l (VI a )(x) \y-x 

JB(xJ){ L j L M-l 

J^B( 



\ 



2M-2 



dy 



V/ ff (y)-P 2M -3(V/ a ,x)(y) 



A u -\VI a ){x) \y-x 



Introduce the function g as above, with P2M-2 replaced by P 2 m-3, so that 



2M-3 



V/ a (y) - Pw-3<yi a ,x)<y) = g(l) - 2 v 



(0) 



f (2M-2)(l-r) 
Jo (2M-2)! 



2M-3 S W ^ 



Now 



(2M-2)! 



I 

^|=2M-2 



c^VUx + r(y - x)) 
)8! 



(y-x) 



= J] (^VI a (x + T(y-x))-d^VI a (x)) 



=2M-2 



(y-xf 



+ 



I 

|=2M-2 



The last term in the preceding equation is the homogeneous Taylor polynomial of de- 
gree 2M - 2 of the vector VI a around the point x. It is then equal to a homogeneous poly- 
nomial of the same degree with zero integral on B(x, t) plus A M ~ l (VI a )(x) \y - x\ 2M ~ 2 
(by the discussion before the statement of Theorem 3 in the introduction). Hence 



J 

JB(x,t) 





ldPVI a {x)\ 


z 

\\B\=2M 


_ 2 l f* J 



)(y-xf--±- 



A M -\VI a )(x) \y-x 



2M-2 



dy = 0, 



and therefore, remarking that f^(2M - 2)(1 - r) 2M 3 dr = 1, 
VK ( t l) (x) 



■ f f (2M-2)(1- 

JB(x,t) JO 



2M-3 



2 (d^VI a (x+T(y-x)) - cfVI a (x)) ^—^ 

V|0|=2Af-2 P' 



drdy . 
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Thus 



|Vtf f (1) (x)| < C f V f teV/aC* + r(v - x)) - ^V/ a (x)) ^— ^ dy 



|/S|=2M-2 

Making the change of variables = r(y - x) the integral in y above becomes 

C hP 
J = T -W \ (cfVI a (x + h)- cfVI a (x)) — dh , 
Jb(0jt) P- 

which is invariant under the change of variables h! = -h, because \/3\ is even. Hence 

hP 



2J = T 



■n 



C hr 
\ (d^VI a (x + h) + d?VI a (x -h)-2 dPVI a (x)) — dh . 

Jb(0jt) P- 



Now we claim that dPVI a is in the Zygmund class for \/3\ = 2M - 2. This follows from the 
fact that the Zygmund class in invariant under homogeneous smooth Calderon -Zygmund 
operators, A M ~' is an elliptic operator and A M ~ l VI a is in the Zygmund class. Hence 



|/| < Cr 



f \h\^ 

JB(0jt) 



dh<Ct 2M - l T. 



Thus 
Since 



(f 

WW/3 



\VK?\x)\ <Ct 2M ~ l 

1/2 



,2(2M-1) 



dt 

+2a+l 



= C \x\ 



-n-l 



we get Hormander's condition with y = 1 in the domain t >\x |/3. 



4.6 The sufficient condition 

In this section we prove that condition (2) in Theorem [3] is sufficient for / e W a ' p . Let 
f,gi, . . .,gN £ L p satisfy S a (f,g\,. . . ,g N ) e LP . Take an infinitely differentiable func- 
tion (p > with compact support in 5(0, 1), J (f> = 1 and set <j) € (x) = pr0(f), e > 0. Con- 
sider the regularized functions f £ = f * (f> £ , gj, £ = gj * (p e , 1 < j < N. We want to show first 
that the infinitely differentiable function f e is in W a ' p . We have (-A) a/2 f e = f * {-Af l2 (p £ . 
We need a lemma. 

Lemma 2. 

(i) If (p is a function in the Schwartz class and a any positive number, then (-A) a/2 (p 
belongs to all U spaces , 1 < q < oo. 

(ii) Iff e L p , 1 < p < oo , then (-A) a/2 f is a tempered distribution. 
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Proof. Set iff = (-A) a/ (p. If a = 2m with m a positive integer, then iff = (—A) m (p is in the 
Schwartz class and so the conclusion in (z) follows. If a = 2m + 1, then 

n 

if/ = (-A) 1/2 (-A)'> = -i £ Rj( dj (-A)"V ) , 

7=1 

where Rj are the Riesz transforms, that is, the Calderon-Zygmund operators whose 
Fourier multiplier is It is clear from the formula above that if/ is infinitely differ- 

entiable on R" and so the integrability issue is only at oo. Since dj (—A) m cp ), has zero 
integral, one has, as x — > oo, \if/(x)\ < C \x\~ n ~\ and so the conclusion follows. 
Assume now that m - 1 < a < m, for some positive integer m. Thus 

If m is even, of the form m = 2M for some positive integer M, then 

^ = A M (p * I m _ a , 

where I,„- a (x) = C \x\ n ' a ~". Hence if/ is infinitely differentiable on R' ! . Since A M <p has zero 
integral, \if/(x)\ < C \x\ n ~ a ~ n ~ l , as x — > oo. But a - m + 1 > and thus if/ is in all U spaces. 
If m is odd, of the form m = 2M + 1 for some non-negative integer M, then 

n 

if/ = -iJ]R j (d j A M i p)*I m ^. 

7=1 

Again iff is infinitely differentiable on R" and, since Rj(djA M (p) has zero integral (just 
look at the Fourier transform and remark that it vanishes at the origin), we get \iff(x)\ < 
C \x\ m ~ a ~ n ~ 1 , as x — > oo, which completes the proof of (i). 

To prove (ii) take a function </? in the Schwartz class. Let q be the exponent conjugate 
to p. Define the action of (-A) a/2 f on the Schwartz function (p as (/, (-A) a,2 (p) . By part 
(z) and Holder's inequality one has 

\((-Af /2 f,<p)\ = \(f,(-AT> 2 <p)\ < C \\f\\ p ||(-Ar/Vll„ 

which completes the proof of (ii). □ 

Let us continue the proof of the sufficiency of condition (2). By the lemma (-A) al2 (f) e 
is in L 1 and so 

iK-Ar /2 / e ii P = ii/ * (-Ar /2 ^ii P < \\f\\p ii(-Ar /2 ^iii . 

Hence f e e W a ' p . 

Next, we claim that 

S a (f e ,gi, e ,...,gN,e)(x)<(S a (f,gi,...,g N )*<f> e )(x), xeR". (23) 
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One has 
where 



S a (f,gl,---, gN)(x) = \\R a (x, t)\\ L 2 Uh/t 2 a+l 



N-l 



R a (x, t) = (f *Xt)(x) - f(x) - Mjgj(x) t 2} - M N (g N *Xt)(x) t 

7=1 



2N 



ifa = 2N and 



N 



R a (x, t) = (f * X t)(x) - f(x) - £ Mjgj(x) & 



7=1 



if 2N < a < 2N + 2 . As before we have set Mj = f \z\ 2j dz . Minkowsky's integral 
inequality now readily yields (1231) . 
Set 

D £ (x) = || 2j Mj (— ^ - 8j , e (x)) - M N |(-^ - g N>e ) * X j (x) t 2N \\ L 2 (dt/t 2^ } 



if a = 2N and 



D £ (x) = || £ Af, - g i>e (*)) ^)|| L 2 Wf2 „ +1) 



7=1 



H2N <a <2N + 2. 
By 



D e (x) < S a (f £ )(x) + SJf e ,gi te , ■ ■ .,g N ,e)(x) 

< S a (f £ )(x) + (S a (f, g!,...,g N )* (p £ ){x) , 

and so D e e V . In particular, D £ (x) is finite for almost all xeR". Thus 



lim inf 

f->0 



N-l 



Mfeix) 



- gjAx)) t 2j - M N I ( — '- - g N , e ) *Xt I (x) t 



7=1 



r ff = o , 



for almost all x e W , if a = 2N , and 



lim inf 

f->0 



7=1 



&f e (x) 

Li 



-gjAx))t 



2j 



r a = o . 



for almost all x € R" , if 2iV < a < 2;V + 2 . It is easy to conclude that the only way this 
may happen is whenever 



A j Mx) 
Lj 



gj, e (x), l<j<N, 
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for almost all x e R". Hence 

A j fe 

Lj 



gj, l<j<N, 



in LP as e — > 0. Since / e — > / in L p as e — > 0, 

A'/ e ->A'/, l<j<N, 
in the weak topology of tempered distributions. Hence 

tJf 



Lj 



gj, l<j<N. 



We claim now that the functions f e are uniformly bounded in W a,p . Indeed, by the 
proof of necessity of condition (2) and by (1231) . 

||(-Ar /2 / e |l P - \\S a (f £ ,Af £ /L u ...,A N f £ /L N )\\ p 

<\\S a (f,Af/L u ...,A N f/L N )\\ p <oo. 

Hence there exist a function h e LP and a sequence e 7 — > as j — > oo such that 

(-Af /2 f £j ^h as j^oo 

in the weak ★ topology of LP. On the other hand, by Lemma 2, {-A) al2 f is a tempered 
distribution and so 

{-Af l2 f e -» (-A) ff/2 / as 6^0 

in the weak topology of tempered distributions. Therefore (-A) a/2 f = h s LP and the 
proof is complete. 

5 Final remarks 

Let (X, d,/j.) be a metric measure space, that is, X is a metric space with distance d and /u 
is a Borel measure on X. We assume that the support of \x is X Then, given a > and 
1 < p < oo, we can define the Sobolev space W a ' p (X) as follows. Let Af be the unique 
integer such that IN < a < IN + 2. Given locally integrable functions f,g\,...,g N define 
a square function by 



S a (f,gi,gz, ...,g N )(.x) 2 = [ f 

Jo JB(xJ) ' 



d/u(y) 



dt 

— , x e 



where D is the diameter of X and i?^(v, x) is 
Rw(y, x) = f(y) - f{x) -gi(x) d(y, x) 2 + gN-i(x) d(y, x) 2(N ~ l) - (gNh( x ,t)d(y, x) 2 
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if a = 2N and 



R N (y, x) = f(y) - fix) - gl (x) d(y, x) 2 + • • • - g N ^(x) d(y, x) 2 ^ - g N (x) d{y, x) 



if IN < a < 2N + 2 . Here the barred integral stands for the mean with respect to on the 
indicated set, B(x, t) is the open ball with center x and radius t and gB(x,t) is the mean of 
the function g on B(x, t). 

We say that a function / belongs to the Sobolev space W" iP (X) provided / e L p {p) 
and there exist functions g u g 2 , . . . , g N e L p Qi) such that S a (f, g u g 2 , ...,g N )e L p (pi). 

We have seen in the previous sections that this definition yields the usual Sobolev 
spaces if X = W is endowed with the Euclidean distance and /i is Lebesgue measure. One 
can prove with some effort that the same is true if W is replaced by a half-space. Very 
likely this should also work for smoothly bounded domains, but we have not gone that 



There are many interesting questions one may ask about these new Sobolev spaces. 
For instance, how do they compare, for a = 1, with the known first order Sobolev spaces, 
notably those introduced by Hajlasz in [0 or the Newtonian spaces of [jSj ? For higher 
orders of smoothness one would like to compare them with those introduced by Liu, Lu 
and Wheeden in ULLWL One may also wonder about their intrinsic properties, namely, 
about versions of the Sobolev imbedding theorem, the Poincare inequality and so on. 

For the Sobolev imbedding theorem the following remark might be useful. In R" 
the LP space can be characterized by means of the following "zero smoothness" square 
function: 



The result is then that a locally integrable function / is in LP if and only if So(f) e LP . The 
proof follows the pattern described several times in this paper. One first deals with the 
case p = 2 via a Fourier transform computation. Then one introduces a L 2 (dt/t) -valued 
operator T such that 



and one shows that its kernel satisfies Hormander's condition. 
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